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The time evolution of O(N) symmetric A$ 4 scalar field theory is studied in the large 
N limit. In this limit the (4>) mean field and two-point correlation function ($$) evolve 
together as a self-consistent closed Hamiltonian system, characterized by a Gaussian density 
matrix. The static part of the effective Hamiltonian defines the True Effective Potential U e f / 
for configurations far from thermal equilibrium. Numerically solving the time evolution 
equations for energy densities corresponding to a quench in the unstable spinodal region, 
we find results quite different from what might be inferred from the equilibrium free energy 
"effective" potential F. Typical time evolutions show effectively irreversible energy flow from 
the coherent mean fields to the quantum fluctuating modes, due to the creation of massless 
Goldstone bosons near threshold. The plasma frequency and collisionlcss damping rate of 
the mean fields are calculated in terms of the particle number density by a linear response 
analysis and compared with the numerical results. Dephasing of the fluctuations leads also 
to the growth of an effective entropy and the transition from quantum to classical behavior 
of the ensemble. In addition to casting some light on fundamental issues of nonequilibrium 
quantum statistical mechanics, the general framework presented in this work may be applied 
to a study of the dynamics of second order phase transitions in a wide variety of Landau- 
Ginsburg systems described by a scalar order parameter. 
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I. INTRODUCTION 

Spontaneous symmetry breaking by a scalar order pa- 
rameter occurs in many different physical systems, as 
diverse as liquid 4 He at temperatures of order 2 K, to 
the standard model of electroweak interactions at tem- 
peratures of order 250 GeV = 3 x I0 15 K. The pro- 
totype renormalizable quantum field theory describing 
this symmetry breaking is a scalar field with a A$ 4 self- 
interaction. The behavior of the finite temperature ef- 
fective potential or Landau-Ginsburg-Helmholtz free en- 
ergy in this theory is well-known and forms the usual 
basis for discussion of the symmetry restoration at high 
temperature. On the other hand, surprisingly little effort 
has been devoted to the nonequilibrium or time depen- 
dent aspects of the symmetry breaking phase transition. 
With the development of practical general techniques for 
studying time dependent problems in quantum field the- 
ory, as well as the advent of high speed supercomputers 
it has become possible to address these dynamical issues 
systematically for the first time. It is clear that a de- 
tailed description of the time dependent dynamics will 
be necessary to calculate nonequilibrium properties of 
the phase transition, such as the formation and evolu- 
tion of defects in the 4 He system after a rapid quench, 
or the efficiency of baryogenesis in the electroweak phase 
transition. Other examples requiring the detailed time 
evolution of scalar fields are the chiral phase transition 
of the strong interactions such as may soon be probed in 
relativistic heavy-ion colliders and the problem of reheat- 



ing the very early universe after it has passed through an 
epoch of rapid expansion and cooling. 

All of these problems require a consistent treatment 
of time-dependent mean fields such as the scalar ex- 
pectation value (3?) in interaction with its own quan- 
tum and/or thermal fluctuations. The technique we rely 
upon in this paper consists of replicating the scalar field 
<fr — > €>i, i — I, . . . , N. Then it becomes possible to com- 
pute the quantum effective action in a systematic power 
series in the parameter 1/N 0. Variation of this effec- 
tive action yields equations of motion for the mean fields 
coupled to Green's functions of the theory which are suit- 
able for implementation on a computer. The leading or- 
der in large N corresponds to a self-consistent mean field 
approximation, i.e. a truncation of the infinite hierar- 
chy of Schwinger-Dyson equations for the n-point corre- 
lation functions to a closed Hamiltonian system of just 
the one-point function, (&(x)) and two-point function, 
(<& (#)<!> (#')). Because no irreducible correlators higher 
than these appear in the leading order of the large N 
expansion, it is equivalent to a Gaussian approximation 
to the time-dependent density matrix of the system. As 
will become apparent, the approximation allows auto- 
matically for a mixed-state Gaussian density matrix, p, 
and is therefore more general than a Gaussian ansatz for 
the pure state wave function in the Schrodinger picture. 
Hence, the mixed-state Gaussian width of p can describe 
at once and on the same footing both the quantum and 
classical statistical fluctuations of the field about its 
mean value far from thermal equilibrium. 
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An important property of the evolution equations in 
the large N limit is that they are Hamilton's equations for 
an effective classical Hamiltonian H e ff (in which Ti ap- 
pears as a parameter). This effective Hamiltonian turns 
out to be nothing else than the expectation value of the 
full quantum Hamiltonian H in the general mixed state 
described by the time-dependent Gaussian density ma- 
trix, p, i.e., 

H eff =Tr(pH). (1.1) 

The canonical variables of H e f f are in one-to-one corre- 
spondence with the parameters needed to specify the gen- 
eral Gaussian density matrix in the Schrodinger picture, 
or the one- and two-point functions of the Schwinger- 
Dyson hierarchy in the Heisenberg picture. In order to 
highlight the Hamiltonian structure of the mean field 
equations, our first purpose in this paper will be to es- 
tablish these various correspondences in detail for the 
O(N) symmetric A<& 4 theory. The existence of an ef- 
fective Hamiltonian makes it clear from the outset that 
the leading order large N approximation is self-consistent 
and energy conserving, and hence does not introduce any 
time irreversibility or dissipation "by hand" into the sys- 
tem. 

A corollary of the identification of the effective Hamil- 
tonian H e ff in the large N limit is that its static piece 
U e ff (obtained by setting all the canonical momenta to 
zero) is the True Effective Potential which governs the 
nonequilibrium evolution of the system. This True Effec- 
tive Potential (TEP) is real, and completely well-defined 
for states far from thermal equilibrium. In the special 
case of thermal equilibrium it becomes the internal en- 
ergy U of the closed Hamiltonian system described by 
H e ff. In contrast, the Hclmholtz free energy F which 
is sometimes called the "effective" potential is not de- 
fined away from thermal equilibrium and becomes com- 
plex in the unstable spinodal co-existence region between 
the two spontaneously broken vacua in simple approxi- 
mation schemes. This is easily understood in the more 
general nonequilibrium context of this paper: it means 
simply that there is no stable equilibrium state in the 
co-existence region with a fixed constant value of (<&). 
In other words, if we start at t = with initial condi- 
tions corresponding to this putative thermal equilibrium 
state, the system immediately begins to evolve in time 
away from it, and this is the physical meaning of the 
imaginary part of F [||. 

We shall see that the free energy function F (or simply 
its real part) is a very poor guide to the time evolution of 
the system far from thermal equilibrium. In particular, 
the oscillations of the time dependent expectation value 
(t) about the spontaneously broken minima are char- 
acterized by a frequency (the plasma frequency) which 
is not the second derivative of F at its minimum (which 
turns out to be zero), and moreover, even the location 
of the minima of F is in general different from the sta- 
tionary points of the mean field evolution. Explicitly 



solving for the actual nonequilibrium mean field dynam- 
ics and demonstrating that it is quite different from what 
might be inferred from an uncritical use of the equilib- 
rium free energy function F is the second major emphasis 
of this work. In this we corroborate the similar conclu- 
sions reached in Ref. ||. 

The spontaneous breakdown of the global O(N) sym- 
metry in the model leads to the existence of TV — 1 mass- 
less Goldstone bosons (in d > 1 spatial dimensions), 
which dominate the dynamics in the large N limit. In 
d < 1 spatial dimension there is no symmetry breaking 
and the system is inevitably driven into the symmetric 
phase, no matter what the initial state or energy den- 
sity. The one dimensional case is of more than passing 
interest in showing how the O(N) symmetry is restored 
dynamically and the Mermin-Wagner-Coleman theorem 
J4[ is satisfied in real time. For d > 1, the absence of 
a finite mass threshold means that an arbitrarily small 
amount of energy in the mean field can create massless 
Goldstone boson pairs nearly at rest. This open chan- 
nel provides an efficient mechanism for the mean field to 
continuously transfer its kinetic energy to the massless 
particle modes over time. The presence of a symmetry 
which requires massless particles is a feature which the 
O(N) model shares with other physically interesting the- 
ories such as non-abelian gauge theories and gravity. The 
O(N) scalar theory provides an instructive example of 
dissipation by means of massless particle creation, which 
should be applicable in other quite diverse contexts, such 
as gluon production in relativistic heavy-ion collisions or 
graviton creation in early universe phase transitions. De- 
veloping techniques and gaining some valuable intuition 
for these more challenging problems is our third reason 
for presenting a study of massless A<E> theory in some 
detail. 

It is remarkable that despite the explicitly Hamiltonian 
structure of the mean field equations, we observe quasi- 
dissipative features in the evolution, in the sense that 
energy flows from the mean field (<&) into the fluctuating 
particle modes without returning over times of physical 
interest. In other words, although the underlying equa- 
tions are fully time-reversal invariant, typical evolutions 
beginning with energy concentrated in the mean fields are 
effectively irreversible, at least over very long intervals of 
time. This apparent irreversibility is quantifiable first, in 
terms of the increase in particle entropy obtained by av- 
eraging over the rapidly varying phases of the fluctuating 
modes, and second, by the effective damping rate of the 
collective motion which we calculate by a standard linear 
response analysis. Since the mean field Gaussian approx- 
imation contains no collision terms, the particles interact- 
ing with each other only through the mean fields, this ef- 
fective relaxation to a quasi-stationary (but non-thermal) 
state is a form of collisionless damping, similar to Lan- 
dau damping in non-relativistic electromagnetic plasmas. 
We call this collisionless damping due to effective loss 
of phase information in the fluctuating quantum modes 
dephasing. Dcphasing of the fluctuations is both an ex- 
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tremely general and efficient mechanism for introducing 
effective dissipation into the reversible Hamiltonian dy- 
namics of mean field evolution. Hence, even in this rel- 
atively simple collisionless approximation to a quantum 
many-body theory one can begin to see how irreversibil- 
ity and the second law of thermodynamics emerge from 
a consistent treatment of fluctuations in a closed Hamil- 
tonian system. Collisions which first appear at one order 
beyond the mean field limit in the large N expansion 
would be expected to make the dephasing and dissipa- 
tion found at lowest order still more efficient. 

In addition to the effective dissipation of energy from 
the collective plasmon mode to the fluctuations, dephas- 
ing is also responsible for quantum decoherence, in the 
sense of suppression with time of the off-diagonal ele- 
ments of the Gaussian density matrix. The point is that 
by going to the appropriate time-dependent number basis 
the diagonal matrix elements of p are adiabatic invariants 
of H e ff and therefore slowly varying functions of time, 
while the off-diagonal elements are very rapidly varying. 
These rapid phase variations in the off-diagonal interfer- 
ence terms cancel out very efficiently when the sum over 
the mode momentum is performed, or if the phases in a 
given mode are averaged in time. In either case, the par- 
ticle creation and interactions have the effect of bringing 
the quantum system into what effectively looks more and 
more nearly like a classical mixture in which the phase 
information in the off-diagonal components may be dis- 
carded for most practical purposes at late times. The 
symmetry breaking behavior of the density matrix and 
effective disappearance of quantum interference between 
the two classically allowed outcomes is the final result. 
An important consequence of decoherence is the appear- 
ance of a diagonal effective density matrix which may 
be sampled to generate smooth classical field configura- 
tions. These configurations are free from spurious cut-off 
dependences and can be used to address issues such as 
the generation of topological defects in a nonequilibrium 
phase transition and the modeling of individual events in 
heavy-ion collisions. 

The time-dependent scalar theory is an excellent theo- 
retical laboratory for the study of general nonequilibrium 
phenomena such as decoherence and the quantum to clas- 
sical transition quite aside from specific potential appli- 
cations to phase transitions in many systems of physical 
interest. The detailed study of effective dissipation and 
decoherence in an explicit field theoretic example is the 
fourth major focus of the present work. Some of our re- 
sults on the Hamiltonian nature of the evolution and on 
dephasing and decoherence have been reported earlier in 
condensed form Q. 

The paper is organized as follows. In the next Section 
we begin by reviewing the large N expansion of A3> 4 the- 
ory in the real time effective action formulation. Then 
the Hamiltonian structure of the equations of motion is 
exhibited and the effective Hamiltonian and Gaussian 
density matrix are identified in Section III. The static 
part of this effective Hamiltonian is the nonequilibrium 



True Effective Potential (TEP), U e ff which we define 
and relate to the thermodynamic free energy F in Sec- 
tion IV. Numerical evolution of the actual equations of 
motion in both one and three dimensions show clearly 
the difference with what might have been inferred from 
F. In Section V we identify the adiabatic particle num- 
ber basis in which dissipation through the increase of the 
effective particle entropy and decoherence are described. 
We also present numerical evidence for the efficient de- 
phasing of the quantum modes in the time-dependent 
mean fields, and show how it leads to typical classical 
configurations in the mixture in which quantum inter- 
ference effects have been washed out. In Section VI we 
perform a linear response analysis of small perturbations 
away from thermal equilibrium, as well as away from the 
non-thermal stationary states found in the previous sec- 
tions, compute the plasmon damping rate, and compare 
it with the numerical results. There is excellent agree- 
ment in the thermal case for both the plasmon frequency 
and damping rate, whereas the non-thermal situation is 
more complicated. We conclude in Section VII with a 
discussion of our results and of their possible application 
to diverse problems of interest in the real time dynam- 
ics of second order phase transitions with a scalar order 
parameter. There are three Appendices, the first on the 
renormalization of the energy and pressure of the the- 
ory, the second cataloguing some mathematical proper- 
ties of the Gaussian density matrix used in the text, and 
the third containing the details of the numerical methods 
used in solving the equations. 

II. THE LARGE N EFFECTIVE ACTION 

The most direct method of deriving the equations of 
motion in the large N approximation is the method of 
the effective action, which also has the advantages of ex- 
hibiting the covariance properties of the theory explicitly 
and providing a general framework for the systematic ex- 
pansion in powers of 1 /N to any desired order beyond the 
mean field approximation. Here we will present a short 
derivation of the effective action and refer the interested 
reader to the earlier work Q for details of the derivation, 
which is quite standard. The underlying O(N) symmet- 
ric scalar field theory with which we begin is described 
in d space dimensions by the classical action, 

Sci{$,x} = / dtd d xL cl [^, X ] 

= J dtd d x |-i$ i G- l [ x ]#i + jx(f +M 2 )} 

(2.1) 

where i — 1 , . . . , N and 

G ~ x \x\ = -n + x, (2.2) 
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with the metric signature (— , +, +, +). This form of the 
action is equivalent to the more familiar Lagrangian den- 
sity, 

1 

,2X 2 



Lci[$] = -z(8v*iW*i)-Voi 



with the definition of the auxiliary field \ by 

A 



X = -At H , 

2N 



(2.3) 



(2.4) 



since the two Lagrangians L c i and L c i are then equal up 
to a surface term. The quartic coupling in the Lagrangian 
has been taken to be A/TV from the outset, rather than 
rescaling it later by 1/N as is sometimes done ||. 

If the parameter /i 2 > 0, then this classical Lagrangian 
describes spontaneous symmetry breaking since the min- 
imum of the classical potential V c i occurs at 



6 iN v vN ; v = \j ^ fi 



(2.5) 



rather than at zero. The second derivative of the classical 
potential at its minimum is 

9 2 Vci A ( ^ x 2N/i 2 \ . A 



d^ l d^ ] 2N 



k^k 



A J N 



0, i or j = 1,...,N- 1 



A«§, i = j = N 



(2.6) 



At this minimum the O(N) symmetry is spontaneously 
broken, x = and there are N —1 massless modes. Small 
oscillations in the remaining i — N direction describe a 
massive mode with bare mass equal to \^2/i = yXvo- In 
this standard way of describing the symmetry breaking, 
one direction is singled out and its vacuum expectation 
value is scaled with \fN as in ( |2.5| ). It is clear that in the 
large N limit the effects of this single degree of freedom 
are down by 1 /N compared to the N — 1 massless modes 
which dominate the dynamics. As we shall see in Section 
VII the effects of the single massive degree of freedom 
are recovered in the behavior of the X correlator at late 
times. 

Passing now to the quantum theory, the large N mean 
field approximation is obtained by retaining only the 
leading terms in a systematic expansion of Feynman di- 
agrams in 1/N. To lowest order, the quantum effective 
action is 

S eff [0, X ] = NS cl [0, x ] + Tr In G" 1 [ X ] . (2.7) 

We use units in which the speed of light and Boltzmann's 
constant are unity, c = fcs = 1, but we retain fi in order 



to exhibit the semiclassical nature of the large N limit. 
The effective action, S e f / generates the proper vertices of 
the 0(N) symmetric A$ 4 theory, correct to leading order 
in 1/N. It is just N times the classical action for a single 
component <$> field, plus the self-consistent (or resummed) 
one-loop correction represented by the last Tr lnG _1 [x] 
term, where x is the mass squared of the propagator G 
appearing in the loops. The quantity <fi = (<&) is the mean 
value of the quantum field 3?. Thus, in this description 
we effectively return to a single component $ theory. 
The role of the large N limit is to justify the neglect 
of the single massive degree of freedom relative to the 
N — 1 massless Goldstone bosons in the propagator G[x] 
of the $ field. However, the massive degree of freedom 
determines the x propagator which controls the plasmon 
oscillations discussed in Section VI. 

Before proceeding it is worth pausing at this poi nt to 
compare and contrast the large N effective action (2.7) 
with two different but related approximations. The sim- 
plest is the standard one-loop approximation on the sin- 
gle component theory, which is obtained by expanding 
the $ field about its mean value <j). The quantum effec- 
tive action in this approximation is 



Si 



-ioo P W\ = Sd[4>] + — TrlnG \V/[] 



where the effective mass appearing in the two-point func- 
tion is 



V", 



d 2 V cl 



a* 2 



3A 



(2.9) 



The simple one-loop approximation includes only the 
one-loop self-energy diagram of Fig. 1, without any fur- 
ther resummation of diagrams. 




FIG. 1. The one-loop self-energy diagram which con- 
tributes to the simple one-loop effective action of Eqn. (2.8) 



In contrast, the approximation usually called the 
Hartree or Gaussian approximation in the literature is 
closely related to the large N approximation in t hat the 
second derivative of the classical potential V c " in ( |2.S| ) is 
replaced by its Gaussian mean value, viz. 



SHartreei^J M 2 ] — S c , 



Tr In G~ l [M 2 } (2.10) 



with 



M 2 



A 

2N 



-V 2 + 7T^7*fe*fe)5ij 



— (***,) 

N 



(2.11) 
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In the Hartree approximation the classical action should 
be expressed in terms of the set of variational parame- 
ters, 4> and Mfy In both the Hartree and large N ap- 
proximations the expectation values of bilincars obey the 
factorization condition, 

(2/)> - (*<(*)} (*i(»)> -iSijG{x,y), (2.12) 

where which two-point function G appears here depends 
on the application. In any case, this factorization implies 
that the last term in ( [2.11 ) above is suppressed by a fac- 



tor oil/N compared to the diagonal 5y term. The dif- 
ference between the Ha rtree and large N effective action 
is just this last term in (2.11). In the large TV approxima- 
tion it is discarded at leading order as in (2.6) and (|2.7|) , 
wher eas in the Hartree approximation of Eqns. ( |2.10| ) 
and (2.11) it is retained. An important consequence of 
this difference between the two approximation schemes 
is that unlike the Hartree approximation, the large N 
effective action is the leading term in a series in a well- 
defined expansion in powers of 1/N. Thus, it is possible 
to improve on the large N limit in a systematic way by 
retaining higher order terms in this series. In contrast, 
the Hartree approximation is simply a variational ansatz. 

In addition to providing a practical expansion tech- 
nique the large N expansion has the significant concep- 
tual advantage over the Hartree approximation of plac- 
ing mean field theory in its proper context with respect 
to other systematic methods in nonequilibrium statisti- 
cal mechanics. In fact, the systematic truncation of the 
Schwinger-Dyson hierarchy of connected 2rt-point func- 
tions by the large N expansion in quantum field theory 
is the precise analog of the truncation of the BBGKY 
hierarchy of n-particle distribution functions in nonequi- 
librium classical statistical mechanics. The existence of 
an energy-conserving expansion parameter in 1 /N which 
preserves all the relevant symmetries of the underlying 
field theory is a powerful technical tool for development 
of the quantum theory of nonequilibrium processes from 
first principles in a systematic way. 

The effective mass which appears in the large N effec- 
tive action is 



= -^ + ^<f> 2 (x)-^G[x}(x,x) 



(2.13) 



which will be recognized as just the expectation value 
of the (operator) definition of the auxi liary field in (2.4) 
upon using the factorization condition ( 2. 12] ). Since G[x] 
itself depends on x through the definition, 



G- 1 \x)°G\ X ] 
(-n+x)G[ X ](x,y) 



1 or 

S d+1 (x- 



V) 



(2.14) 



in position space, ( 2.1 3| ) is a nonlinear integral equation 
for the auxiliary field x- The Feynman diagrams that 
contribute in both the large N limit and the Hartree ap- 
proximation are the sum over all daisy and superdaisy 



diagrams represented in Fig. 2, with the only difference 
that the lines correspond to the propagator G[x] in one 
case (large N) and G[M 2 ] in the other case (Hartree). 
This is just the resummation of self-energy diagrams re- 
quired by a renormalization group analysis, and in fact, 
X is a renormalization group invariant. In the large N 
limit the fundamental excitations are the N— 1 Goldstone 
modes, whose masslessness is fixed by the 0(N) symme- 
try which the 1 /N expansion respects order by order. In 
contrast, since the Hartree approximation is not a sys- 
tematic expansion in any small parameter it is not clear 
a priori how to renormalize M 2 , nor whether it has any 
renormalization group invariant meaning. The result is 
that although the two approximations are very similar 
in some respects and may be handled by the same tech- 
niques, the physics they describe is really quite different. 
Which approximation is more reliable depends very much 
on the application, and in particular whether or not the 
massless Goldstone modes of the large N limit actually 
play the lead role in the physics we wish to describe. 




FIG. 2. Typical daisy and superdaisy diagrams included 
in the Hartree or leading order large N approximation to the 
effective action. 

In the real time formulation, the effective action of 
( |2.7| ) becomes the starting point for all further analysis 
of th e large N dynamics. From the construction of iS in 
( |2.7D as a Legendre transform in the presence of external 
sources, it follows that its first variation with respect to 
the independent mean fields ^ or x is proportional to 
the sources for these fields. In the absence of external 
sources, these variations yield the equations of motion 
for the mean fields, 



and 

1 6S e ff[(f>,x] _ (x(x) 



N 



Sx 



(2.15) 



= o 



(2.16) 



which just recovers ( 2.13 ) to this order. Unlike the equa- 
tion ( [2.15 ) for (j>, ( 2.16 J involves no derivatives and is 
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therefore an equation of constraint (or gap equation), 
rather than an equation of motion for an independent 
propagating degree of freedom. That x 1S nevertheless 
a useful indicator of symmetry breaking shoul d be c lear 
from the spacetime independent form of Eqn. ( [2.15 ), 



(d„4> = 0) 



(2.17) 



which tells us that cither cj> or x ( or both) must vanish in 
a uniform, stationary state. The case, <j> = is the O(N) 
symmetric state with generally positive Xi while the case 
of non-vanishing <fi is the spontaneously broken state in 
which x = is the vanishing Goldstone boson mass. This 
is an explicit realization of Goldstone's theorem, which 
is respected by the large N approximation. 

In the following the explicit functional dependence of 
G[x] and its inverse on the mean field x will be sup- 
pressed, and we adopt the simpler notation G(x, y) or 
still more briefly G hereafter. 

The important difference of the large N e quatio ns from 
the purely classical ones is the last term in ( 2.13 ) involv- 
ing TiG. Whereas the spontaneous symmetry breaking 
solution x — can be achieved only when <f> — ±vq = 
±y^2/Xfj, classically, in the large N approximation the 
two-point function G contributes at the same order and 
can even dominate the mean field <f> in Eqn. (2.13). Since 



G itself depends on x this additional term also has the 
effect that the x field can undergo nonlinear collective 
plasmon-like oscillations as we shall see explicitly in Sec- 
tion VII. Notice also that there is nothing to prevent x 
from being negative at some times which allows us to 
explore the dynamics of the unstable spinodal region of 
Fig. 3. The spinodal at a given temperature is the region 
where the second derivative of the potential becomes neg- 
ative. In the large N approximation, this corresponds to 
X<0. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

*/v 

FIG. 3. The spinodal value of <f>, </>t is the value of <j> 
where x nrs t becomes negative. It satisfies the equation 
(<t> T /v) 2 = 1 - (T/T c ) 2 shown above (see Eqn. (O) below). 



In this paper we shall be concerned only with spatially 
homogeneous mean fields <fi = <fi(t) and x — xW- It is not 
difficult to treat spatially inhomogeneous mean fields by 
the same methods, but as the homogeneous case is sim- 
pler and already contains much of the essential physics, 



we restrict ourselves to that case in this paper. When <f> 
and x are functions only of t, then the two-point Green's 
function G is a function only of the spatial difference 
x — x' and it is useful to introduce the Fourier transform 
of the corresponding Wightman function G> (or G<), 

G > (f,x;i',x') = -G*(t,x;t',x / ) 

= y[dk]e lk -( x - x ')G>(M';k) (2.18) 



with 



[dk] 



d d k 

{2n) d 



(2.19) 



In fact, G> is a function only of k = |k| by rota- 
tional symmetry of the spatially homogeneous state. The 
Wightman functions are solutions of the source- free wave 
equation G _1 o G> = G _1 o G< = 0. In terms of complex 
mode functions fk(t) obeying 



A(t) = o, 



(2.20) 



it is convenient to express the Wightman functions as 

G>(i,i';fc) = if k (t)f* k (t')(N(k) + l) + if* k (t)fk(t')N(k) , 

(2.21) 

where N(k) (not to be confused with the N of large N) 
is an arbitrary time-independent function of k. It car- 
ries the interpretation of particle number density in the 
general spatially homogeneous initial state in the basis 
specified by the mode functions and is arbitrary ex- 
cept for the requirement of being fi nitely integrable with 
respect to the integration measure (2.19). 

Explicitly, if the original quantum field $ is expanded 
about its mean value in terms of the quantum modes 
and the corresponding plane wave creation and destruc- 
tion operators 4 and ak, 



*(t,x) = (*> +£*(*, x) 



= 0(t)+L-^(e ik - x / fc (i)a k 



in a cubical box of finite length L, then 
N(k = |k|) = (4 flk ) 



(2.22 



(2.23) 



is just the expectation value of the particle number den- 
sity in this basis. It is also a function only of k by isotropy 
of the spatially homogeneous state. The Wightman func- 
tions are given then by the usual expressions, 

hG> (t, x; t', x') = i{6Q(t, x)5$(t', x')) , 
aG<(i,x;t',x') = i(»(t',x')<5$(t,x)) , (2.24) 



G 



and the condition G> = — G< follows immediately from 
t he h crmiticity of the underlying quantum field (5<& in 
(2.22). The commutation relation 



<5k,k' 



(2.25) 



implies that the complex mode functions should be cho- 
sen to satisfy the Wronskian condition, 



f $1 _ f*^!± 
Jk dt Jk dt 



ih . 



(2.26) 



in order for the quantum field operator <fr to obey the 
usual canonical commutation relation, 

([#(*, x), ^(*,x')]> = ^^(x-x') . (2.27) 

The norm alizat ion of the time-independent Wronskian 
condition ( 2.26| ) is the only place where the constant h 
enters the mean field equations, which otherwise are quite 
classical in their time evolution dynamics. 

We have not writte n the a priori possible bilinear terms 
fkfk or /fc/jj? in ( 2.21 ), since they can always be absorbed 
into a redefinition of fk and N(k) under the transforma- 
tion, 

fk - cosh 7fe e ie " +i ^f k + sinh 7fc e **-<**/jJ (2.28) 



without a ffectin g the mode equation (2.20) or Wronskian 
condition ( 2.26 ). This is equivalent to making a Bogoli- 
ubov transformation which sets to zero the expectation 
values of the pair densities (a-^a-^) = (a^a^) = 0. Hence 
there is a natural SU(l, 1) Heisenberg group structure 
(for each k) inherent in the leading order large N equa- 
tions. 

In the gap equation ( 2.1 3| ) the coincidence limit of the 
Green's function G(x, x) appears. The coincidence limit 
of either G> , G< or the Feynman propagator G are all 
identical, so we obtain from Eqns. (2.21) and (2.18) the 



coincidence limit of any of the Green's functions in the 
form of an integral over k of the Fourier mode functions, 

fiGy (t, x; t, x) = fiG < (t, x; t, x) 

= i f [dk]\f k (t)\ 2 (2JV(A) + 1). (2.29) 



Since the mode functions obeying ( [2.20 ) and (2.26) be- 
have as 



Mt)^ yC|Ie X p(-f/fc(tO) (2-30) 



for large k, where 

w fc (t)= (k 2 + x(t)Y 



(2.31) 



the integral in ( 2.29| ) is quadratically divergent in d = 3 
spatial dimensions. Introducing an explicit momentum 



cut-off A and perfo rming the angular k integrations, the 
gap equation ( 2.15 ) may be rewritten in the form, 



k 3 dk |/ t (()| 2 c„ 



for d = 3, where we have defined the notation, 
a k = 2N(k) + 1 . 



(2.32) 



(2.33) 



The fact that the bare parameters and Aa must de- 
pend on the cut-off in order to render the equations in- 
dependent of A in the end has been exhibited explicitly 
as well. The quadratic divergence in d = 3 is the di- 
vergence of the one-loop self energy diagram in Fig. 1 
and is absorbed into the bare mass parameter fi\. One 
convenien t way to effect this mass renormalization is to 
evaluate ( 2.32| ) in the time-independent spontaneously 
broken vacuum, \ — 0: 4> — v i thereby absorbing the 
quadratic divergence into the relation between the bare 
and physical expectation value of the field. In this way 
we eliminate ji\ and obtain 



x(t) 



Aa 
'Air 2 



■A / k " dk \ \Mt)\ 2 ^ k --\, (2.34) 



which is free of quadratic divergences. 

The re maini ng logarithmic divergence in the mode in- 
tegral of (2.34) is removed by the logarithmic coupling 
constant renormalization in the usual way, i.e. 



A A = ^ 1 (A,m) X R (m 2 ) 



(2.35) 



with 



Z x (A,m) = 1 



Ai?(m ) In — 



32tt 2 

h 



32tt : 



rA A In 



A 2 



(2.36) 



and \n(m 2 ) the renormalized <fr 4 coupling defined at 
some finite mass scale to 2 . By dividing both sides of 
( [2.34 ) by Aa and using Eqns. (2.35) and (2.36), it is 



straightforward to verify from the large k behavior of the 
integrand t hat t he logarithmic dependence on A of th e 
integral in ( |2.34| ) is cancelled by the logarithm in ( 2.36 ). 
Thus, the resulting equation for %(i) is independent of 
A for A large, and x is m fact a renormalization group 
invariant physical mass squared of the theory. The con- 
dition that Z\ > prevents us from taking the cut-off 
strictly to infinity with Xr > fixed, for otherwise Z\ 
from ( 2.36] ) would eventually become negative and the 
theory would become unstable. This is just a reflection 
of the Landau ghost instability of scalar $ field theory, 
and means that the theory can be sensible and non-trivial 
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only as an effective field theory equipped with a large but 
finite cut-off, A. This presents no problem in practice as 
long as A is large enough that the physical time evolu- 
tion, plasma oscillations, damping, etc. occur on time 
scales much greater than A -1 . In that case, the evolu- 
tion is numerically quite insensitive to the value of A over 
a very wide range, provided A_r is not too large Q|. 

There is no wave function or <p renormalization at low- 
est order in large N so no further renorma lizat ion is re - 
quire d in d = 3 to this order, and Eqns. ( 2. IE ), ( 2.20 ), 
and ( |2.34D together with the constraint on the initial 
data Q2.26| ) specify a well-defined closed system of evolu- 
tion equations for the mean fields <p and \ m interaction 
with the fluctuations Let us emphasize again that 
these equations differ from the purely classical tree level 
approximation or the simp le o ne-lo op approximation in 
that the last term of ( 2.32| ) or ( [2.34 ) couples the fluctua- 
tions self-consistently and nonlinearly back on the time- 
dependent mass gap function x(£), corresponding to the 
full sum of daisy and superdaisy diagrams in Fi g. 2. 

In d = 1 space dimensions the integral in (2.13) is 
only logarithmically ultraviolet divergent and (2.34) is 



already well-defined and finite, without any A renormal- 
ization. However, in lower spatial dimensions the small 
k or infrared behavior of the integral becomes more deli- 
cate, and must be treated carefully. 



III. THE EFFECTIVE HAMILTONIAN AND 
DENSITY MATRIX 

The presentation of the large N equations of mo- 
tion of the previous section was based on the functional 
method, in which the extremization of the effective ac- 
tion S e ff in ( 2.1 5| ) and ( 2.1 6| ) takes the place of the usual 
Euler-Lagrange variational principle of classical mechan- 
ics. The equations so obtained are for the mean values of 
the field operators and their two-point Green's functions 
which describe fluctuations about the mean fields. This 
immediately raises a question: Is there a corresponding 
Hamilton form of the variational principle for an effec- 
tive large N Hamiltonian involving both the mean fields 
and their fluctuations? In addition one would like to 
know: To what distribution of field amplitudes in the 
Schrodinger wave function (or density matrix) do the 
large N equations for the one and two-point functions 
correspond? It is to these questions that we turn in this 
section. Answering them will lead directly to the true 
effective potential of nonequilibrium large N mean field 
theory. 

Let us begin by consideration of the case of d = 
spatial dimensions, i.e. quantum mechanics. The gener- 
alization to higher d will turn out to be s trai ghtfo rwa rd. 
For d — the Lagrangian L c / or L c i of (2T) or (2J3) is 
that of an A" component anharmonic oscillator. If A" = 1 
it reduces to the usual anharmonic double well oscillator. 
In the d — case there is no k index to be integrated and 
the equations derived in the last section become simply 



dt 2 



+ X (t) ) cb(t) = 
A 



X(t)= g (4> 2 (t)+e(t)-v 2 ), (d = 0) (3.1) 

where we have introduced the notation, 

e(t)^a\f(t)\ 2 ^(2N+l)\f(t)\ 2 (3.2) 

in terms of the expectation of the number operator, A" = 
(a^a). We recognize that 

e{t) = mm)) w = (m) - m 2 ) m 

is just the quadratic variance of the quantum variable 
<&(£) from its mean value (f>(t). By differentiating this 
relation and defining rj = £ we find 



mut) = ((*(*)*(*) 

= 2a\Re(//*) . 



•#(t)*(t) - 2</>(t)j>(t) 



(3.4) 



The advantage of defining these quantities will become 
apparent from the role they play in the physical interpre- 
tation of the equations of motion (|3.l|) i n the Schrodinger 
picture. Indeed, by differentiating ^. 4|) again and using 
both the equation of motion (2.20) (with k = in d = 
spatial dimensions) and the Wronskian condition (2.2£) 
for f(t), we find 



rj = £ 



2„2 



-xV 



l£3 



(3.5) 



in the ne w notation. This equation of motion together 
with (3.1) are just Hamilton's equations, 



dH, 



P = 



V 



eff 



dH, 



eff 



8^ 



(3.6) 



for the effective iwo-dimensional classical Hamiltonian, 

+^ 2 +e-vir + ^f (3.7) 



with 



P = <t> 
rj = i- 



dH, 



eff 



dp 
dH eff 
drj 



(3.8) 



the canonical momenta conjugate to the two generalized 
coordinates 4> and £. A different but equivalent set of 
canonical variables was discussed in Ref . ||7| . 



Hence by the simple change of notation in (3.2) and 



(3.4) we have recognized that the large A^ equations for 
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the quantum anharmonic oscillator with classical poten- 
tial 



A 

8iV 



JV 



^ I ! 



(3.9) 



derived by the effective action technique of the last sec- 
tion are precisely equivalen t to Hamilton's equations for 
the effective Hamiltonian (3.7). This answers the first 
question we posed at the beginning of this section in the 
affirmative, at least for the case of d = 0. 

An immediate corollary of the Hamiltonian structure 
in the extended phase space (p, </>;7?,£) is that the large 
N evolution equations are energy conserving with H e ff 
the value of the conserved energy. It is also interesting 
to note that ha is a constant of the motion which enters 
the effective potential 



U, 



eff\ 



2^2 



(3.10) 



together as a single "centrifugal barrier" term, whose ef- 
fect is to repel the variance £ away from zero. The phys- 
ical and mathematical analogy to an angular momen- 
tum barrier is made even stronger by the fact that the 
three symmetric bilinears aa, aW and aa^ + a^a gen- 
erate the Lie algebra of su(l, 1) or so(2, 1) which is the 
non-compact version of the ordinary angular momentum 
algebra, su(2) or so(3), and moreover, the Casimir in- 
variant of this rank one Lie algebra is exactly ?i 2 (T 2 /4, 
in the standard normalization. The corresponding Lie 
group is just the three parameter group of homogeneous 
linear Bogoliubov transformations ( [2.2S ) . 

To answer the second question posed at the beginning 
of this section let us recall that both the time-dependent 
Hartree (TDH) and large N equations have been studied 
in the Schrodinger representation, and they are known to 
correspond to a Gaussian trial wave function ansatz Q . 
Indeed, it is straightforward to verify that the Gaussian 
ansatz for the normalized pure state Schrodinger wave 
function, 



(27r^(i))-3exp<^ 



.p(t)x 



1 , ■ '/" = \ ,,,,2 



4£ 2 (t) 2H(t)J 



{x-m\ (3.ii) 



obeys the expectation value of the Schrodinger equation, 



<*(*)! -tE 



n 2 " a 2 



2 ^ d&id®i 

2 — 1 

d 



vm i*(*)> 



(3.12) 



in the coordinate representation where, with P as the 
canonical momentum, 



(*(t)|0(*,P)|*(t)) 

= J dx$*(x;t)0 (x,-iJi^J *(x;t) , (3.13) 

provided t he l arge N lim it is taken and the equations 
of motion (3.1) and (|3.5|) are satisfied for a = 1. Thus, 



the Gaussian ansatz (3.11) is a special case of the general 
large N equations of motion, where £ and r\ are related to 
the real and imaginary parts respectively of the Gaussian 
covariance. 

In earlier work Q it had been recognized that the 
Gaussian ansatz for the Schrodinger wave function(al) 
imposed one constraint on the ab initio three indepen- 
dent symmetrized variances, 

(*(t)|(*-0) 2 |*(t)) 
(*(i)|(P* + *P 2#)|*(i)) 



2£7? 



(*(t)|(P-p) 2 |*(t)> = f? 2 



4£ 2 



(3.14) 
(pure state), 



expressing all three in terms of only the two variables £ 
and 77, in the present notation. The one anti-symmetrized 
variance is fixed by the commutation relation, [#,P] = 
ih. To what does the restriction to a = 1 correspond 
and how can it be relaxed? The answer to this question 
i s sug gest ed b y the form of ( [3.14 ), and definition of a in 
( |2.33 ) or (3.2), which shows that a = 1 corresponds to 
zero expectation of the number operator a^a, i.e. to the 
pure state vacuum annihilated by a. However, the mean 
field equations of the Section II allow for the more general 
possibility that this expectation value may take on any 
constant value N. For example, we might consider the 
finite temperature Bose-Einstein distribution, 



iV T = 



exp 



a, = I +2V, =coth ( ^ ] I 



(3.15) 



for some loq and temperature T. Since such a thermal 
state corresponds not to a pure state Schrodinger wave 
function, but rather to a mixed state density matrix, it 
is not surprising that a pure state Gaussian wave func- 
tion ansatz cannot describe this case. However, the gen- 
eral structure of the mean field equations involves only 
the one and two-point functions of the quantum variable, 
so we should expect them to correspond to a Gaussian 
ansatz but for a mixed-state density matrix instead of a 
pure state wave function. The most general form for this 
mixed-state normalized Gaussian density matrix p is 

(x'\p(p,<l);r]^;a)\x) = (2tt£ 2 )-* expj i |(a;' - x) 



a + 1 \i ' 



8£ 2 

JL 

'2h£_ 



[(x' - 0) 2 + (x - 



<7 Z -1 



(3.16) 
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in the coordinate representation. In the special case that 
(7 = 1 the last (mixed) term in the exponent vanishes and 
p reduces to the pure state product, 



in the large N limit where e is the energy density defined 
in ( |A2| ). The three symmetrized variances are indeed now 
all independent with 



P (t) =mt))m)\ 



(3.17) 



with \^(t)} given by ( |3.11 ). For a > 1 the general Gaus- 
sian p does not decompose into a product, and 



/oo />oo 
dx I dx'(x\p(t)\x')(x'\p(t)\a 
-OO J — OO 



= a- 1 <l 



(3.18) 



which is char acteristic of a mixed state density matrix. 

That ( fc.l6|) is indeed the correct gene ralization of the 
pure state Gaussian wave function (3.11) is easily verified 
by checking that p(t) satisfies the expectation value of the 
quantum Liouvillc equation, 



Tr [ift^p) = Tr [H. /i] 



(3.19) 



provided the large N limit as befo re is taken and the 
equations of motion (3.1) and (3.5) are satisfied for ar- 
bitrary a. Taking the large N limit is equivalent here to 
the replacement of the full anharmonic Hamiltonian by 
a time-dependent harmonic oscillator Hamiltonian, 



H -> H 



,,^i(P 2 +c 2 (t)* 2 ) 



(3.20) 



where 



d 2 V 



\^ 2 {t)+e{t)-vi)= X {t) 

(3.21) 



is the self-consistently determined frequency of the oscil- 
lator in the large N limit. With this replacement it is 
straightforward to verify that the Gaussian form is pre- 
served by the time evolution under H osc (l(]] ]ll]] ]l2] |, 



(3.22) 



In fact, subsitution of the Gaussian form (3.16) into this 
Liouville equation and equating coefficients of x, x', x 2 , 
x and xx' gives five evolution equations for the five 
parameters speci fyin g the Gau ssian which are none other 
than the Eqns. (3.1) and ( |3 . 5| ) together with a — 0. 

The effective classical Hamiltonian for the large N 
equations, ( |3.7P is just the expectation value of the quan- 
tum Hamiltonian in this general mixed state Gaussian 
density matrix, i.e. 

H eff (<j>,p;Z,n;a) = Tr (pH) = Tr ( P H osc ) = eL d , 

(3.23) 



Tr(p(*-0) 2 ) =e , 
Tr (p (P* + *P 2<jyp)) = 2^ , 

Tr(p(P-p) 2 ) = V 2 



h 2 a 2 
4^2 



(3.24) 



replacing (3.14) of the pure state case. The mean values, 



<j> = (*) = Tr (*p) ; and 



(*) = Tr (Pp) 
(3.25) 



remain valid for both the pure and mixed state cases. 

The physical interpretation of the five parameters of 
the general large N equations (p,<fi,r),!;;a) in terms of 
the general time-dependent mixed-state Gaussian den- 
sity matrix of the Schrodinger picture is now explicit in 
d = quantum mechanics. Since a is the constant pa- 
rameter which determines the degree of mixing and h 
and a appear only in the combination ha it is clear that 
the large N equations allow for a smooth interpolation 
between the quantum pure state case in which ha = h to 
the high temperature or classical limit where 



2T 

har — > — 

UJQ 



T->oo or h^O , (3.26) 



in which h drops out entirely. Thus, quantum and clas- 
sical thermal fluctuations are treated on the same foot- 
ing in the large N limit, with the value of the constant 
parameter ha determining whether the fluctuations de- 
scribed by £ and ry are to be regarded as predominantly 
quantal or thermal, or intermediate between the two. 
Moreover, we see that large N Gaussian dynamics is re- 
ally classical dynamics of a Gaussian distribution func- 
tion, except that ha which measures the second moment 
of the classical distribution cannot be taken to zero as it 
could be classically, but instead is bounded from below 
by h. This is made explicit by the form of t he W igner 
function corresponding to the density matrix ( 3. 16 ) , viz. 



fw{x,p x ) 
1 



1 



irha 



exp 



2nh _ 
I 2£ 2 



dye 



-ip x y/h 



■ l y I I 

x+^\p\x- 



2£ 2 



1 2 



px-<i>- j{x- 



(3.27) 



Before leaving our d — example it is instructive to 
examine the static solutions of the effective Hamiltonian, 
viz. the simultaneous vanishing of 



P 



dH, 



eff 



X 



b 2 +i 2 -vl)<j> = x^ = ® , and 



h 2 - 2 



4£ 3 



= 



(3.28) 
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If we look for a spontaneosly broken solution <j> ^ 0, then 
X must vanish from the first of these conditions. But then 
we cannot satisfy the second condition for finite ha and 
£. This is just a rederivation of the fact that there can 
be no spontaneous symmetry breaking in d — quantum 
mechanics (or in fact for any d < 1). We are forced 
instead to the symmetry restored situation for which <f> = 
and x = Xo > is determined from the real positive 
root of the cubic equation 



Xo = ^(£o ~ v a) 



2 ^2 



(3.29) 



for Co( fT )- The Gaussian density matrix centered at 
= with variance £, 2 (<j) is the solution of the time- 
independent Liouville equation for the anharmonic dou- 
ble well oscillator with a trial variational density ma- 
trix of the form (3.16). In the limit that the height of 
the energy barrier between the two wells, Et, — Awq/8 is 
much greater than the fluctuation energy in either well, 
Ef = hay/Xvo/2, the width of the Gaussian £g — > Vq, 
and its energy E -> h 2 a 2 /8v% = Ej/l6E b < E f , corre- 
sponding to a probability density spread over the entire 
region (— vo, vq) in all N components of <E>i. 

The entire development of the Hamiltonian equations 
and Gaussian density matrix is quite easy to general- 
ize to any number of spatial dimensions d, at least for 
the case of spatially homogeneous mean fields. Since in 
Fourier space the mode equations are just replicated at 
every spatial momentum k we have simply to introduce 
the subscript k on all of the relevant definitions in this 
section. For example, the density matrix for d > can 
be written as a product of Gaussians in Fourier space, 
viz. 

({<PkY\p\Ww}) = Y[(Wk}'\p(Pk,(l)k;Vk,S,w;crw))\Ww}) 

k 

= Il( 2 <k)^exp{ l ^(^-^ k ) 

k ^ 

" 1 [Ok - M 2 + (vk - k ) 2 ] 



2^k 



4£ 



-Wk - 0k)(<^k - ^k) 



(3.30) 



where tp^ is the generalized coordinate of the field am- 
plitude in Fourier space and pk the corresponding canon- 
ical momentum. The mean fields and their canonical 
momenta 



0k = S k0 <j)(t), p k = $kop(t) 



(3.31) 



all vanish except for k = in the spatially homogeneous 
case. The definitions 



(3.32) 



and 



x(*) = ^(0 2 (*) + 



<f(t)+ [dk] urn 



h 

2k 



(3.33) 



have been introduced in obvious analogy to the d = 
case. The effective Hamiltonian density which gives rise 
to these equations is 



H, 



eff 



L d 



1 2 
£ = -p 



1 

2A ; 



[dk] 



Hi 

(3.34) 



with x regarded as a depend ent v ariable of and the 
£fc through the gap equation ( 3.33 ) above. The k 2 term 
arising from the spatial gradient of the field in d > 
dimensions is the most significant difference from the d = 
case considered previously. The Hamiltonian equations 
of motion 



P , P 



dH, 



eff 



£,k =Vk , Vk = - - 



dH, 



eff 



96 



-X0 , 

-{k 2 + x)ik + n2al 



(3.35) 



are completely equiva lent t o the equations of motio n for 
t he m ode functions ( 2.20 ) and mean fields Jill ) and 
( |2.13 ) in the spatially homogeneous case in any number 
of dimensions. 

This completes our demonstration that the large N 
mean field equations of the previous section are Hamil- 
ton's equations with H e ff given explicitly by (3.34) 
above, and the identifi catio n of the time-dependent Gaus- 
sian density matrix ( 3.30 ) to which the homogeneous 
mean field equations correspond. The case of spatially 
inhomogeneous mean fields involves only the straight- 
forward generalization of (3.30) to Gaussian covariances 
which arc off-diagonal in the momentum index k, and a 
corresponding coupling between the diffe rent momentum 
modes in the effective Hamiltonian (3.34). Since we have 
no need of these expressions in the present work we do 
not give the explicit formulae here, though they may be 
easily worked out within the present framework. 



IV. THE NONEQUILIBRIUM TRUE EFFECTIVE 
POTENTIAL (TEP) 



Having obtained the effective Hamiltonian which de- 
scribes the evolution of the closed system of mean fields 
and fluctuations in the large N limit, it is natural to de- 
fine the dynamical or true effective potential to be the 
static part of the effective Hamiltonian density (3.34), 
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Ueff (<!>,{&} ;W) 



2A' 

X 
2 



83 



(4.1) 



S 2 - v 2 - - 



\ I [dk] (V 



4$ 



This [7 e / / is the energy d ensity of an initial state with a 
Gaussian density matrix (3.3C) centered around cj> with 
instantaneously zero velocities p — r)k — in both the 
mean field and the fluctuation variables. We call U e ff 
the "True Effective Potential" (TEP) because it deter- 
mines the true out-of-equilibrium time evolution dynam- 
ics of the system according to Hamilton's Eqns. (3.35). 



It must be clearly distinguished from the finite tempera- 
ture free energy effective potential often discussed in the 
literature. The relationship between the two we would 
like to discuss next. 

The standard method of calculating the free energy at 
temperature T is to analytically continue the effective ac- 
tion S e ff to imaginary time |l3| and evaluate the Tr In 
G _1 term over fluctuations which are periodic in imag- 
inary time with period /3 = K/T. In this way one finds 
for the case at hand, 



i(* 



Tin 



1 — exp 



T 



which is a real function of <f> provided 
u\ = k 2 + X > . 



(4.2) 



(4.3) 



The infinite zero point energy in F is handled in the usual 
way by subtracting the energy of the zero temperature 
vacuum at x = 0, i.e. f[d]i]hk/2. With the subtrac- 
tion of this temperature independent constant and the 
previous definitions of v and the logarithmic coupling 
renormalization by Eqns. (2.35) and (2.36), the ther- 



modynamic free energy becomes finite and is given by 



F«.r) = ^(* 2 -» 2 -^f) 



-2 Tin 



1 — exp 



flU>k 



(4.4) 



with x T the solution of the finite temperature gap equa- 
tion, 

ftA f 00 „ f i , fhuj k \ 1} , . 



As long as x T > the free energy function (4.4) is real 
and well-defined. 



by 



The first derivative of the free energy function is given 



dF 



(4.6) 



upon using the gap equation (4.5). In the spontaneously 
broken state defined by x T = 0, we have 



v 2 - 



h 

2^ jo 
12% ' 



k dk 



exp (M) - i 



(4.7) 



The last relation informs us that the expectation value 
4>t vanishes at the critical temperature, 



T c = 2V3hv. 



(4.8) 



At T — T c , the second derivative of F((f>,T) vanishes at 
= 0, and therefore the phase transition is second or- 
der in the large N approximation. We remark that this is 
not the case in the simple one-loop or Hartree approxima- 
tion where the transition is weakly first order [[15| . Since 
the <E> 4 field theory lies in the same universality class as 
spin models which are known to have second order phase 
transitions, the large N mean field approximation gets 
the correct order of the transition, and is another reason 
why it is to be preferred over the other approximations 
Jig] . Indeed, the correlation length, which is given by the 
the inverse mass of the radial excitation, 



t(T) = \ R *4>j 



( 6 



V xrT c 



(4.9) 



diverges as T — > T c with the — i critical exponent of 
mean field theory. 

Actually the second derivative of the free energy func- 
tion vanishes at its minimum for any T <T C . T his fo llows 
the fact that the derivative of the gap equation (4.5) eval- 
uated at x T — involves an infrared linear divergence so 
that 



dx T {<t>) 



0. 



x T =o 



and 



(4.10) 



(4.11) 



at x T =0. This perhaps surprising feature of the large 
N free energy is illustrated in Fig. 4. It is a direct 
consequence of the massless Goldstone particles present 
in the spontaneously broken phase. As we shall see in 
Section VII, it also makes the equilibrium free energy 
useless for describing the plasma oscillations about the 
spontaneously broken thermal minimum. 
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FIG. 4. The Helmholtz Free Energy F as a function of (f> 
for three temperatures T = .1, 1, 1.5 (in descending order 
from the top) in units of T c in the leading order large N ap- 
proximation. The curves terminate at the minimum of F at 
4> — (fir, XT = 0, given by Eqn. (4.7), which is the boundary 
of the spinodal region, i.e. for cf> < cj>T there is no stable equi- 
librium state. The second derivative of F at <f> — 4>T vanishes, 
showing in particular that the phase transition at T = T c is 
second order. 



For T < T c and 4> < 4>t the gap equation (4.5) has 
no positive real solutions. This is the unstable spinodal 
region of Fig. 3. If % < then the frequency u>k becomes 
imaginary for k < |x| 3 - This means that strictly speak- 
ing the thermodynamic free energy F is not well defined 
in this region. Despite the fact that no thermal equi- 
librium uniform state can exist in this spinodal region, 
what is sometimes done is to define the free energy by an 
analytic continuation from X > to x = ~lxl < 0. By 
this procedure F acquires an imaginary part, while its 
real part is no longer a convex function of <j) (as general 
theorems require) This has led to much discussion in the 
literature |17j |18|. However, we needn't wonder at the 
meaning of this analytic continuation and imaginary part 
of F in the more general nonequilibrium context. The an- 
alytic continuation from x > to x = — Ixl < simply 
defines a Gaussian density matrix centered at <f> < 4>t 
with instantaneously zero velocities, <j> — rjk = and 
parameters and Ok — 2Nx(k) + 1 determined by the 
thermal distribution at \\ T \ — ~X T where the free energy 
is well defined. The density matrix will certainly evolve 
away from this unstable configuration and the imaginary 
part is a measure of the rate at which this evolution will 
occur initially However, shortly afterwards the non- 
linear effects of the evolution described by our general 
nonequilibrium equations of motion will set in, and this 
imaginary part can only be at best a simple order of mag- 
nitude estimate for the evolution away from the initial 
state in the weak coupling limit \r <C 1. 

The nonconvexity of the real part of the free energy is 
also related to the same instability |l8j . The point is that 
the construction of the equilibrium free energy by a loop 
expansion or the large N expansion tacitly assumes the 



existence of a stable configuration as the starting point 
for the expansion. It is precisely this assumption which 
breaks down in the unstable spinodal region. A careful 
definition of the "effective potential" F as the minimum 
of the free energy for fixed <fr = ($) independent of space 
and time leads to a flat constant function between the 
two minima at ±v which is convex in accordance with 
the general theorems. However this minimization cannot 
be ac hieved with a Gaussian density matrix of the form 
Q3.30D and the flat convex form of F defined in this man- 
ner tells us nothing about the dynamical evolution of an 
initial Gaussian centered at a value of cj> in the spinodal 
region. For this reason we propose to drop all attempts 
to refine the definition of a free energy "effective poten- 
tial" and focus instead on the potential which actually 
governs the nonequilibrium dynamics in a given approxi- 
mation scheme. In the leading order large N scheme this 
dynamical potential is the True Effective Potential U e f / . 

In contrast to the free energy F, the TEP U e ff is a 
function(al) of all of the generalized fluctuation coordi- 
nates £fc as well as the mean field 0. It also depends on 
the constant parameters o^, w hich need not be that in 
the thermal distribution (3.15). Hence it is defined much 
more generally than F and is manifestly real and posi- 
tive. If we wish to consider a function only of the mean 
field <p it is possible to try first minimizing U e ft with 
respect to the Gaussian parameters 



au, 



eff 



2 ^2 



= 



which has the real solution 



p2 _ TlVk 



(4.12) 



(4.13) 



provided again that (4.3) holds 



If we restrict U e ff further by requiring Ok to be the 
Bose- Einstein thermal distribution, <JTk — l + 2JVr(fe) of 
Q3.15Q we find 



u{<t>,T) = E/e//0M6fc} = {6J ;CT fc = <7 Tk ) 

X ( ,2 2 X N 

' [dk]Lo k {l + 2N T (k)) . (4.14) 

The first term is the classical potential energy density V c i 
of the mean field <p and the second term is the quantum 
plus t hermal energy of the fluctuati ons a t the extremum 
(4.13). We recognize the expression (4.14) as the in terna l 
ener gy of the Gaussian configuration specified by ( 4.13| ) 
and (3.15), which differs from the Helmholtz free energy 
F by the standard thermodynamic relation, 



F=U -TS 



(4.15) 



with 
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S = — Tr pt InpT 

= J [dk] [(N T (k) + 1) In (JV T (fc) + 1) - N T (k) lnN T (k)] 

(4.16) 

the von Neumann entropy of a Gaussian thermal density 
matrix pT- 
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FIG. 5. The True Effective Potential U e ff as a function of 
<j> and evaluated in a thermal equilibrium state for the same 
three temperatures as in the previous figure (note here that 
the lowest temperature is now the bottom curve). The Gaus- 
sian width para mete r has been set equal to its minimum 
value by Eqn. (4.12). Although qualitatively similar to the 
free energy F in Fig. 4, U e ff becomes equal to the internal 
energy U in thermal equilibrium. 

Although F and U have qualitatively similar behav- 
iors in the region \ > where they are both defined 
unambiguously, as shown in Figs. 4 and 5, clearly F and 
U have quite different physical meanings and applicabil- 
ity. The first is the negative of the pressure of the gas 
of scalar particles against which work must be done in 
compressing the system at fixed temperature, i.e. if the 
compression is performed while in contact with an exter- 
nal heat reservoir at temperature T. The second is the 
energy density of the closed system which is conserved if 
it is allowed to evolve, isolated from all external sources 
or sinks of energy. The free energy F can be defined only 
in thermal equilibrium. The internal energy U is a spe- 
cial case of the more general true effective potential U e f / 
which is defined for any Gaussian density matrix, equi- 
librium or not, and which determines the evolution of 
the system away from the instantaneous stationary state 
where < j> = tjk = according to the Hamiltonian equa- 
tions, J3.35 ) . Moreover, in the spinodal region where 
T <T C and (f> < <j>T it is clear that the minimization con- 



dition (4.12) cannot be satisfied for real £?, which informs 



us that the time derivative of the canonical momentum, 
f]k > for k 2 < \x\, i.e. the configuration is unsta- 
ble against the fluctuation widths growing in time for 
these low k. Note that by the definition ( 3.32| ), in terms 
of the mode function is necessarily real. The phys- 



ical interpretation is thus quite straightforward in the 
more general nonequilibrium framework we have laid out: 
there simply is no stationary Gaussian density matrix for 
X < 0, and any such initial configuration will necessarily 
evolve in time by growth of the small k Fourier com- 
ponents of the Gaussian parameters Indeed, this is 
obvious from the mode equations (2.20) which show that 
all the long wavelength modes with k l < \\\ will grow ex- 
ponentially rather than oscillate when \ < j|. Notice 
that the non-existence of a solution to the minimization 
condition ( 4.1 2| ) implies that the exponentially growing 
instability lies in the fluctuations for small k. Only 
the nonlinear back reaction of this exponential g rowth 
of the modes in the time-dependent gap equation ( 2.13j ) 
can eventually bring \ to non-negative values and turn 
off the instability. It is clear that this time-dependent 
process involves many low k Fourier components of the 
field <& and cannot be described adequately by a single 
function of one variable such as the "effective" potential 
F, much less a function of one variable defined only in 
thermal equilibrium in the spinodal region where no such 
thermal equilibrium exists. 

That the time evolution as determined by U e t / is quite 
different from what might be inferred from the free en- 
ergy function F is illustrated by our numerical solution 
of the evolution equations, presented in Fig. 6. The 
oscillation frequency and even the final point of the evo- 
lution of cj> as t — ► oo in general bear no simple relation 
to the minimum <pT of the free energy "effective" poten- 
tial. This is easily understood from the stationary points 
of the TEP, for a true stationary state does exist in the 
spontaneously broken phase if we also require 



dU, 



eff 



d4> 



X0 = O 



since this is satisfied by 



x = o 



and 



ft 



~2k 



(4.17) 



(4.18) 



Notice that there is such a static so lutio n for any <f> and 
<7fc. The finite temperature case, (3.15) is only one of 
infinitely many possibilities for a static solution of the 
mean field equations. Correspondingly there are also an 
infinite number of stationary density matrices satisfying 
[p, Hose] = 0, one for each choice of a^. If Ofc = l + 2JV(fc) 
is given, then <j) is determined by the static condition 
X = 0, i.e. 



lx=o 



h / [dk] 



N(k) 



h 

2^2 



kdkN{k). 
(4.19) 



For arbitrary N(k) the relation ( 4.19 ) may be viewed 
as a kind of sum rule which allows us to distribute any 
fraction of the v 2 in the coherent vacuum expectation 
value (f> 2 and the remainder in the integral over particle 
fluctuation modes. In particular, one may even consider 
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the case where the mean field <f> = identically, but the 
theory is still in the spontane ously broken phase since 
X vanishes and the integral in ( 4.1 9| ) saturates v 2 . This 
case is a kind of marginal symmetry breaking since any 
small perturbation of 4> will take us away from zero mean 
field. 



With the help of the sum rule (4.19) one can easily 
understand the perhaps initially surprising result, first 
found numerically in Ref . , that the evolution of the 
(f> mean field can settle to a value different from, and 
even much smaller than the thermal equilibrium value 
4>t- It simply corresponds to the fact that the distribu- 
tion of particles in the final state is not at all a thermal 
one. Obviously such a situation cannot be described by 
the thermodynamic free energy F since its evaluation by 
continuation to imaginary time assumes from the out- 
set a thermal distribution of particles. For this reason 
one could not expect the true evolution according to the 
mean field equations of motion to bear much relation to 
what one might infer from an uncritical use of F, par- 
ticularly in the region <f> < <j)x where no thermodynamic 
uniform equilibrium state exists and the Hclmholtz free 
energy is not even strictly defined. 
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FIG. 6. A typical evolution of the mean field (j> and the 
mass squared \ in the spontaneously broken phase, starting 
from an unstable initial state with x(0) — —0.37 < 0. All 
time and lengths are measured in units of v = 1. Note that 
if we use the initial energy density e = 0.0697 to infer an 
effective temperature of massless bosons, then the expected 
minimum of F occurs at 4>out = 0.98 which is considerably 
different than the observed <f> ut = 0.81. The oscillation fre- 
quency about <f> — (p ou t also bears no relation to the second 
derivative of F which vanishes at its minimum. 

The inclusion of collisions in the next order of the ex- 
pansion in 1/N makes it more likely to bring the dis- 



tribution of particles closer to a thermal one. However, 
since the evolution is unitary, this can happen only in 
some effective sense. This aspect is closely connected to 
the notion of dephasing which will be discussed in detail 
in Section VI. 

Finally we point out that the stationary spontaneously 
broken solution \ — is disallowed in d = 1 spatial di- 
mension (as it was in d = quantum mechanics), since 
the integration in the definition of x m ( |3.33| ) diverges 
logarithmically as k — > which is inconsistent with x 
vanishing in (4. IS). It is this infrared divergence which 
prevents spontaneous symmetry breaking of the global 
O(N) symmetry in one space dimension, consistent with 
the Mermin-Wagner-Coleman theorem Q. In two or 
higher space dimensions there is no such divergenc e and 
the spontaneous symmetry breaking static solution (4. IS) 
is allowed. 
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FIG. 7. Evolution of the mean field <j> and the mass 
squared \ m one space dimension, starting from an unsta- 
ble initial state. The middle figure shows that \ does not go 
to zero asymptotically, and the oscillations do not damp as 
effectively as in three dimensions. Both these behaviors are 
consequences of the Mermin-Wagner-Coleman theorem which 
prevents spontaneous symmetry breaking in one space dimen- 



The consequences of the Mermin-Wagner-Coleman 
theorem for noncquilibrium dynamics in one dimension 
are illustrated in Fig. 7 where the numerical evolution of 
X and <p is shown beginning from an unstable initial state 
with negative x- We see that x does not go to zero at 
late times (no symmetry breaking) and that both the x 
and 4> oscillations damp very slowly as compared to the 
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behavior in three dimensions. 



V. DEPHASING, DISSIPATION AND 
DECOHERENCE 

Although we have shown by explicit construction of 
the effective Hamiltonian (3.34) that the large N mean 
field equations are completely time reversible, neverthe- 
less one observes in typical evolutions as in Fig. 6 an 
effective irreversibility, in the sense that energy flows 
from the mean fields 4> and x to the fluctuating modes 
fk without returning over times of physical interest. It 
is our purpose in this section to explain this apparent 
irreversibility in terms of dephasing, i.e. the dynamical 
averaging to zero in the sums over k of the rapidly varying 
phases of the fluctuations at a given time. To the extent 
that this phase averaging is exact and the information 
in the phases cannot be recovered, the time evolution is 
irreversible. This is the sense in which the fluctuations 
fk act as a "heat bath" or "environment" for the mean 
field evolution of <f> and \- Of course, since informa- 
tion is never truly "lost" in a closed Hamiltonian system 
evolved with arbitrarily high accuracy, the information in 
the phases can be recovered in principle and we should 
expect Poincare recurrences in the mean field evolution 
after very long times. As the number of modes fk (and 
particularly as the number of relevant infrared modes) 
approaches infinity, we would expect the recurrence time 
to go to infinity. In typical evolutions we followed several 
tens of thousands of modes, and recurrences were never 
observed in practice. The precise dependence of the re- 
currence time on the number of modes is an interesting 
question which may be studied quantitatively within our 
mean field framework by numerical methods. However, 
we have not undertaken such a systematic study here, 
and leave the question of the recurrence time for future 
research. In any case, we cannot expect the leading or- 
der large N collisionless approximation to continue to be 
valid for times longer than the collisional relaxation time 
in the full theory. The study of recurrence times in the 
leading order approximation may be interesting neverthe- 
less as a model for how such recurrences and rephasing 
can occur in more realistic situations. 

To understand precisely what is meant by dephasing 
let us return to the Gau ssian density matrix of only one 
degree of freedom ( 3.16 ). The coordinate representation 
is only one of many possible representations of the den- 
sity matrix. The number basis, defined by the integer 
eigenvalues of 



a' a\n) = n 



(5.1) 



defines a time- independent basis associated with the Fock 
decomposition of the Heisenberg operators (for d = 0), 



$(t) = ( f>(t)+af(t) + a^r(t), 
*(t) = P(t)=^(t)+o/(t)+oV*(*) 



in terms of its mean value and qu antize d fluctuations. 
By using the Wronskian condition ( 2.26| ) we may solve 
these relations for a and a', 



(*-</>)/*- (P - 0)/* =iha, 
(* - <P)f - (P - 4>)f = -iha) . 



(5.3) 



In this time independent basis, the Hamiltonian is given 
by 



H = hu(t = 0) [a)a+\ 



(5.4) 



This transformation is just a (complex) canonical trans- 
formation which on the quantum level is implemented 
by a unitary transformation of bases, with the transfor- 
mation matrix, (x\n) in Dirac's notation. To find this 
transformation matrix explicitly we use the definitions 
( |3.2[ ) and (3.4), the equation of motion (3.5), and the 
Wronskian condition again to secure 



h 2 a ( n\ 



1 



aa+ 2 



c 2 (p-p) 2 



£,((#-0(p- P ) + (p-p)(#-0) 



in terms of £, ?y, and a . Then by acting with this operator 
identity on the matrix (x\n) and replacing — > x and 
P — * —ifa-j^. we obtain the differential equation, 

Ti 2 <j(n + —)(x\n) = 



(5.5) 



which is easily solved in terms of the ordinary harmonic 
oscillator wavefunctions, 




by 

(x\n) = exp 



1 



(2 n n\)i W?i 



—x | exp 



ip(x — <p) irj{x — <f>y 



2h£ 



(5.6) 
(5.7) 



With the transformation matrix element determined 
it is a straightforward exercise in Gaussian integration 
and the properties of Her mite p olynomials to obtain the 
form of the density matrix ( 3.16| ) in the time- independent 
number representation, 



(n'\p\n) 



dx 



dx' (n | x ) (x \p | x) (x | n) 



2<L 



(5.2) 



(5.8) 
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where a is defined by (2.33). The derivation of these 
results and some further properties on the transformation 
matrices are given in Appendix B. 

The matrix elements of the Liouville equation in this 
basis are 



in which 



ih{n\p\ri) = (n\[H,p]\n') 

= {E n - E n ,)(n\p\n') 



(5.9) 



Clearly, if the density matrix in this basis is initially diag- 
onal, which can always be accomplished by a Bogoliubov 
transformation (Appendix B), then it stays diagonal and 
time-independent. The time independence of the matrix 
elements of p in the Heisenberg number basis \n) is simply 
a reflection of the fact that the density matrix (like the 
state vector |^;)) is time independent in the Heisenberg 
picture where th e fiel d operator $ depends on time ac- 
cording to Eqn. (2.22). All the time-dependent dynamics 
resides in the transition matrix element (x\n) which de- 
pends on the five variables (4>(t),p(t); £{t),ri(t); a) while 
the matrix elements (5.8) remain forever unchanging un- 
der the mean field evolution. Indeed, this is just another 
reflection of the unitary Hamiltonian nature of the evo- 
lution, since the von Neumann entropy of the Gaussian 
density matrix 



H osc = — {a'a + aa 1 ) 



(5.15) 



is diagonal for each time. In the a' a number basis, p is 
no longer diagonal, (a), (act), etc., are non- vanishing, and 
A*" = (a^a) 7^ A in general, except in the static case of 
constant uj. 

In addition to diagonalizing the time-dependent har- 
monic oscillator Hamiltonian H osc which describes the 
evolution of the Gaussian density matrix, the h basis 
has another important property, namely the existence of 
an adiabatic invariant W . This adiabatic invariant may 
be constructed from the Hamilton- Jacobi equation cor- 
responding to the effective classical Hamiltonian H e f f in 
the usual way, 



1 (dWX 

2 Kite J 



i fdwy 

2 



2 ^2 



E. 



(5.16) 



This equation is separable and the ansatz W(x, £) = 
Wi(x) + W 2 (0 along with E = E x + E 2 yields, in the 
first variable, 



S = — Tr p In p 



i 



a- 1 



2 J V 2 j V 2 
(A + l)ln(A + 1) - N In N 



In 



a - 1 



(5.10) 



is strictly constant under the time evolution. Hence, 
there is no information lost in the evolution in any strict 
sense. 

There is however another basis which is more appro- 
priate for discussing "physical" particle number. This is 
the time-dependent Fock basis specified by replacing the 
mode function / which satisfies 



by the adiabatic mode function, 



— exp (-i ^ dMt>) ) , (5.12) 



/(*) 



and the Fock representation ( |5 .2[ ) by 

#(*) = a(t)/(t) + at(t)/*(t), 

*(t) = P(t) = -iuj(t)a(t)f(t) + iuj(t)^{t)f*(t) , (5.13) 

where the a and operators must now be time- 
dependent. The corresponding time-dependent number 
basis is defined by 



a^a|n) = n |n) , 



(5.14) 



Wi = jdx V2^T 



while in the second variable ^ one obtains 



W 2 



dt X 2E 2 -^e 



h z a 2 
l£2~ 



(5.17) 



(5.18) 



Under the substitution r — (, 2 /2 this turns out to be 
the same integral which occurs in the Kepler problem, 
again pointing out the formal similarity with the angular 
momentum barrier in a central potential. Using standard 
methods lEOf, one finds 



Wo = 2tt- 



E 9 



Hence the full adiabatic invariant is 
W 



Iixn nio 2 



N . 



(5.19) 



(5.20) 



upon using the definition of A, the definition of a and 



E 



(H) 



(H„ 



/,._,> ! N + ~ 



(5.21) 



Since A is strictly a constant of the motion, Eqn. ( 5.20| ) 
informs us that A is an adiabatic invariant and there- 
fore slowly varying if lo 2 = x is a slowly varying function 
of time. In the language of classical action-angle vari- 
ables, W is an action variable which is slowly varying 
while the angle variable conjugate to it varies rapidly in 
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time (linearly with t ime i n the limit lu is a constant, as 
in the exponent in ( 5.12 )). This means that although 
the density matrix p is certainly not diagonal in the h 
number basis, its off-diagonal elements which depend on 
the angle variable will be very rapidly varying functions 
of time, whereas its diagonal elements in this basis will 
be only slowly varying. This is clearly seen in the ex- 
plicit form for the matrix elements of this basis in Eqn. 
(B22) of Appendix B. Hence if we are interested only in 
the motion of the mean fields which are slowly varying 
functions of time we may average over the rapid phase 
variations in the off-diagonal matrix elements of p in this 
adiabatic number basis. For only a single quantum de- 
gree of freedom this amounts to a time-averaging, and 
can be implemented only by fiat. However in field theory 
there are many momentum modes, so that this averaging 
is actually performed for us in the mean field evolution 
equations by the integrations over k at fixed t. 




















FIG. 8. Integrand and integral of the RHS of Eqn. ( |5.23| ) 
as a function of k for fixed t = 257. 
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FIG. 9. Integrand and integral of the phase averaged quan- 
tity defined by the adiabatic particle number basis as a func- 
tion of k for the same fixed t = 257. The value of the integral 
using the replacement (5.22) is the same as that of the previ- 
ous figure. 

The term dephasing has a precise meaning in this adi- 



abatic particle number representation of the density ma- 
trix. To the extent that the phases in the off-diagonal ma- 
trix elements of the density matrix in the adiabatic h ba- 
sis are rapidly varying in time, they should have little or 
no effect on the evolution of more slowly varying quanti- 
ties such as the mean fields 4> and \ or t ne mean adiabatic 
particle number (N) itself. Thus, a natural approxima- 
tion to the full density matrix p is immediately suggested 
by the existence of the adiabatic numbe r bas is, namely 
to discard the off-diagonal elements (see (|B22j ) below) of 
p in this basis, which is equivalent to time-averaging the 
distribution of fluctuations in the exact Gaussian density 
matrix (3.3C) over times long compared to their rapid 
variations. In field theory, where there are many Fourier 
modes, each with its own rapid phase variation, the ef- 
fect of dephasing may be obtained by simply integrating 
over the momentum index k at fixed time. In either case, 
we expect this averaging procedure to scarcely affect the 
actual evolution of the mean fields, and the extent to 
which this expectation is realized is the extent to which 
dephasing of the fluctuations is effective, and the evolu- 
tion is irreversible. In Figs. 8 and 9 we compare the 
actual k dependence of the function G(k, t) appearing in 
the mean field evolution equation ( p. 13 ) with the phase 
averaged quantity defined by the adiabatic particle num- 
ber N(k), i.e. we make the replacement, 

-ihG(k, t) = \f k (t)\ 2 -> ^(2N(k) + 1) . (5.22) 

The r elevant integral in the renormalized gap equation 
(HI is 



G{x, x, x) - G(x, x,x = 0) = 

5,jf ^{lAWI^-^}. 



(5.23) 



Notice that although there are many mor e oscillations 
in the integrand of the RHS of Eqn. (5.23) than in the 
phase averaged quantity, the integrals over k of the two 
quantities are indistinguishable. 

The adiabatic particle basis may also be used to de- 
fine effective Boltzmann and von Neum ann e ntropies pl| ] 
through the diagonal matrix elements ( B22 ) below. Ne- 
glecting the rapidly varying off-diagonal matrix elements 
gives 



Sb = 



{ (ft $) + i) ln (ft $) + l) - ft(k) lniV(fc)} 



Seff = -Tr p eff ln p eff 



^2^2p2e{k) In p 2 i(k) 
k i=o 



(5.24) 



for the Boltzmann and effective von Neumann entropy 
of the truncated density matrix, where p2i(k) is given by 
Eqn. ( 5.25| ) below. The evolution of these two quantities 
for a typical solution of the mean field equations is shown 
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in Fig. 10. Both display general increase due to continu- 
ous creation of massless Goldstone particles near thresh- 
old 1 22 1 . Neither quantity is a strictly monotonic function 
of time and neither obeys a strict Boltzmann ii-theorem. 
Since the particle modes fk interact with the mean field 
X but not directly with each other, the effective damp- 
ing observed is certainly collisionless, and the dephasing 
here is similar to that responsible for Landau damping 
of collective modes in classical electromagnetic plasmas. 
The entropy S e f / of the effective density matrix provides 
a precise measure of the information lost by treating the 
phases as random. The Boltzmann "entropy" would be 
expected to equal S e f / only in true thermodynamic equi- 
librium, which is not a chieved in the collisionless approxi- 
mation of Eqns. ( 2.20 ). Notice the non-thermal distribu- 
tion of particles in Fig. 9. In this non-thermal distribu- 
tion we see from Fig. 10 that the Boltzmann entropy Sb 
generally overestimates the amount of information lost 
by phase averaging. 
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FIG. 10. Evolution of the Boltzmann and effective von 
Neumann entropies of the diagonal elements of the density 
matrix in the adiabatic particle number basis. 

To the extent that the phase information in the off- 
diagonal matrix elements of p is irretrievable the system 
has become effectively classical, in the sense that the 
quantum interference effects present in the original en- 
semble represented by p are washed out as well. In that 
case we might as well regard the ensemble represented 
by the diagonal, truncated effective density matrix in the 
adiabatic number basis as a classical probability distri- 
bution with the diagonal elements of (h = 2£\p\n = 21) 
giving the classical probabilities of observing £ pairs in 
the ensemble. This probability distribution is derived in 
Appendix B and given by 



P2i(k) = {n % = 2£\p\n % = 21) 
{2£-l)\\ 



~ 2 i e] : sech 7 ^ tanh 2f 7 fe 
where 7 is the magnitude of the Bogoliubov transforma- 



(5.25) 



tion between the a and a bases, given explicitly in terms 
of the mode functions by 



A^) = sinh 2 7? = lfk + lUJkfkl 



2HuJk 



(5.26) 



which depends only on the magnitude of k by spatial 
homogeneity. Sampling this distribution with random 
phases will yield typical classical field amplitudes which 
make up the distribution. Obtaining such typical classi- 
cal fields in the ensemble can give us explicit realizations 
of the symmetry breaking behavior of the system, as well 
as providing the starting point for the study of topologi- 
cal defects produced during the phase transition. 
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FIG. 11. Four typical field configurations drawn from the 
same classical distribution of probabilitie s in the adiabatic 
particle number basis, according to Eqn. (3.30), for t=490 in 
the case the mean field cj> = 0. The units of r are u _1 and 
Aa = 1. 

The sampling of the field configurations proceeds in 
several steps. For a fixed late time and Fourier wave 
number k we calculate the Bogoliubov transformation 
coefficient from the static n to time-dependent adiabatic 
particle number n basis from Eqn. (5.26). This gives a 
set of numbers p^i = (n = 2£\p\h — 2£) normalized to 
unit total probability, 



(5.27) 



and which typically fall off very rapidly with £ so that 
only a finite number of the pit need be retained. Then 
we sample this distribution by drawing a random number 
q in the unit interval [0,1]. Looking at the table of pn 
and the partial sums, 



(5.28) 



e'=o 



we find £ such that 
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Qt-i <q<Qe (5.29) 

to determine h = 2£ of this random drawing for the given 
value of k. We then write 

^ *) = -w J2 ~w= ( flkeik ' x + a k^ lk,x ) 

— > — ^- [ dk k sin krJ — — cos 6 k , (5.30) 
n 2 r J V 2cj fc 

where we have performed the angular integrations in d = 
3 dimensions, and written 

flfc = >/^e <9 * (5.31) 

in terms of a random pha se After performing the 
Fourier transform in ( ]5.30D the result is a typical field 
configuration as a function of radial r at the fixed time 
t. In this way we obtained the four field configurations 
shown in Fig. 11. We observe that although the mean 
field (j> = when averaged over the entire ensemble, typi- 
cal field configurations in the ensemble are quite far from 
zero. In fact they sample values of <fi between the two 
minima at v and —v. We observe as well that there is a 
typical correlation length in the classical field configura- 
tions which is of the order of the inverse momentum in 
which the power of the two-point function is distributed, 
as in Fig. 9. 

It is clearly possible by such sampling techniques to 
generate typical field configurations with any number of 
components $^ in any number of spatial dimensions d. 
Then by appropriately matching the number of compo- 
nents to the number of dimensions we could search for 
different types of defect-like structures such as vortices, 
strings, and domain walls. This classical description of 
the quantum mean field theory is possible only when 
a classical decoherent limit exists through the diagonal 
density matrix in the adiabatic particle number basis. 
Moreover, the classical fields generated by this proce- 
dure are smooth and free of any cut-off dependent short 
distance effects. The extraction of smooth classical field 
configurations from a quantum mean field description is 
a direct consequence of dephasing and the definition of 
defect number in these smooth configurations is free of 
the difficulties encountered when classical definitions of 
defects are applied uncritically to quantum field theories. 
It would be very interesting to pursue these ideas further 
in more realistic models of phase transitions where such 
topological defects are expected. This we leave for a fu- 
ture investigation. 

The dephasing of the density matrix justifies the re- 
placement of the exact Gaussian p by its diagonal ele- 
ments only in the adiabatic number basis leading to the 
effective density matrix p e f f . Now if this effective density 
matrix is transformed back into the coordinate basis it 
corresponds to a density matrix of the original Gaussian 
form but with zero mean fields, zero momentum rjk = 
and <7fc = l + 2A(fc) replacing a, i.e. we obtain a product 
of Gaussians of the form, 



{<Pk\Peff\<Pk) = 

(2«£)-i expj-g^ - <p k f - ^{fp' k + ^ fe ) 2 | , 

(5.32) 

in all the k ^ modes where the mean field vanishes. 
This form shows that the Gaussian distribution off the 
diagonal ip' k = (p k is strongly suppressed compared to the 
diagonal distribution. Indeed in the numerical evolutions 
we have found that kN(k) behaves like a constant for 
small k. This means that a cx 1/k which is large for 
small fc. Hence whereas the root mean square deviation 

from zero along the diagonal ip' k = ip k is = \J T%N(k)/k 
which is large, the root mean square deviation from zero 
in the orthogonal direction off the diagonal is 

^= / h 
<Jk V 2k N(k) 

which is much smaller. The result is that there is vir- 
tually no support in this distribution for "Schrodinger 
cat" states in which quantum interference effects between 
the two classically allowed macroscopic states at v and 
—v can be observed. Instead we see in a different way 
how dephasing produces an ensemble which may be re- 
garded as a classical probability distribution over classi- 
cally distinct outcomes (i.e. the diagonal field amplitude 
distribution) but with essentially no components in the 
off-diagonal quantum interference between these classi- 
cal configurations. This is illustrated in Fig. 12. Hence 
even in this rather simple collisionlcss approximation the 
particle creation effects in the time dependent mean field 
give rise to strong suppression of quantum interference 
effects and mediate the quantum to classical transition 
of the ensemble. 




FIG. 12. The Gaussian p e ff for k = .4 for the data and 
parameter values of Fig. 9 illustrating the strong suppression 
of off-diagonal components due to dephasing. 

When more realistic collisional interactions are in- 
cluded at the next and higher orders beyond the leading 
mean field limit we can expect this transition to decoher- 
ent classical behavior to be even more pronounced. As 
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the Gaussian assumption is relaxed we expect the single 
peak at the origin to split into two peaks. Hence we can 
begin to see how a quantum phase transition leads to an 
effective classically broken symmetry in which large do- 
mains are in a definite classical ground state or another 
but not in a quantum superposition of ground states. 



VI. LINEAR RESPONSE, PLASMA 
OSCILLATIONS AND DAMPING RATE 

The very efficient damping of the oscillations around 
the final state which we have studied in the density ma- 
trix formalism of the last section may also be understood 
as due to the continuous creation of massless Goldstonc 
particles near threshold. In this appearance of strictly 
massless modes the spontaneously broken A$ 4 model is 
qualitatively different from our previous studies of QED 
in the large N approximation |23|. In that case the 
charged particles are massive and there is finally a tunnel- 
ing barrier which shuts off particle creation effects after 
the mean electric field has decreased below a certain crit- 
ical value. Beyond this point the mean field undergoes 
essentially undamped plasma oscillations since there is 
no transfer of energy to created particles. In contrast, 
the present model has no such critical threshold in the 
mean field x which is free to continue creating massless 
bosons for arbitrarily small amplitude and \ asymp- 
totically. This very efficient asymptotic damping of the 
X mean field towards its stationary spontaneously broken 
solution at x = is very well illustrated by the numerical 
results shown in Figs. 13 and 14 for Aa = 1 and an initial 
X = -.5 




0.1 - 

X o - 

-0.1 - 
-0.2 
-0.3 

-0.4 - 



FIG. 13. Time evolution and effective damping of the \ 
mean field towards zero, its stationary value in the sponta- 
neously broken phase. 

Since x becomes very small at late times, it is possi- 
ble to analyze the approach to zero by linear response 
methods. In principle this is possible for any initial con- 
dition but we restrict ourselves to the thermal and vac- 
uum cases. With thermal initial conditions, both the fre- 



quency and the damping rate of the plasma oscillations 
can be calculated and are in very good agreement with 
the simulations as will be discussed further below. The 
long time behavior for initial conditions with N(k) = 
is somewhat more complicated. While the oscillation fre- 
quency is still very well determined by the linear response 
analysis, the oscillation envelope does not damp exponen- 
tially. Instead it falls off as an exponential multiplied by 
a power law in time. 




FIG. 14. Same evolution as in Fig. 13 but followed to large 
times and with an expanded scale for \ to show the long time 
behavior of the time-dependent Goldstone mass squared. 

Given any static solution to the equations of m otion 
with spontaneous symmetry breaking, as in ( 1.1§| ), we 
may consider the real time linear response of the system 
away from this static solution. This is accomplished by 
linearizing the evolution equations in the deviations from 
the static solution, 



x o + x 

fk -» fk(t) + 6f k (t) 



with 



fk(t) = \l — exp(-ikt) 



(6.1) 



(6.2) 



the mode functions corresponding to the static solution, 
X = 0. 

The linearized mode equation, 



= -X(t)fk(t) (6.3) 

may be solved by use of the retarded Green's function, 

„ , , sin kit — t') , ,. , 
G R (t-t';k) = K - >-6(t-t') (6.4) 

in the form, 



dt 2 



Sfk(t) = A k f k (t) + B k fi(t) 



dt'G R (t-t';k) X (t')f k (t') (6.5) 
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where Ak and Bk are coefficients of the solutions to the 
homo gene ous equation. Because the Wronskian condi- 
tion (2.26) is maintained under the linearization, the Ak 
must have vanishing real part, 

ReA k = 0. (6.6) 

The linearized <f> equation gives in the same way, 

8<t>(t) = t 8<j>(0) + 8(f>(0) 

dt' G R {t -t';k = 0) X (t')(f>(f) • (6-7) 

At the same time the linearized gap equation reads 

X = \^8cf> + \ f[dk]Re(8f k f k *)(l + 2N(k)) . (6.8) 



Upon substituting Eqns. (6~5|-3.7) this becomes a linear 
integral equation for the perturbation, x, 



x(t) = -X f dt' n(t - t') X {t') + XB(t) 
Jo 



(6.9) 



where 

IL(t) = t(j) 2 



[dk] ^ + 2 , y k " sin U cos fc* (6.10) 



is the polarization part in the static background and 
B(t) = t(f) 8(f)(0) + (f>8(f>(0) 

[dk] (2+^m Re{Bke2lkt) (611) 

depends only on the initial perturbation away from the 
static solution. The linearized integral equation ( |6.9| ) 
may be put in the form 



dt' D- l (t-t') X (t') = -B(t) 



(6.12) 



where D^ 1 the x inverse propagator function. 

The most direct method of solving such integral equa- 
tions is to make use of the Laplace transforms, 



fl(s) ee 



dte~ st 



U(t) 

(l + 2N(k)) 
2k(s 2 + 4fc 2 ) 



B(s) 



X(s) 



dte' st B(t) 



cf, 8(f>(0) 



m <I±f2» Re 



B, 



s — 2ik 



dte~ st x(t) = -D(s)B(s) , 



(6.13) 



and D(s), the simple algebraic reciprocal of the Laplace 
transform of D _1 (t), 



~ 1 , s 1 =w s 1 <f> 2 
-D-\s) = - x+ Il( S ) = - + ^ 



[dk] 



(l + 2JV(fc)) 
2fc(s 2 + 4fc 2 ) ■ 

(6.14) 



Both the bare coupling A = Aa and the integral over k 
involving the N(k) independent term are logarithmically 
divergent. They combine to give the ultraviolet finite 
contribution, 



1 



1 



Afl(s 2 /4) A A 32^ 2 
1 



h , /4A 2 



h ( 4m 2 
In 



\ R (m 2 ) 32tt 2 V s 



(6.15) 



The renormalized sum rule ( fOj| ) may be used in the 
N (k) dependent term to secure as well 



1 



Afl(5 2 /4) 



h 

2^ 



kdkN(k) 



(6.16) 



which is now independent of the static value of < 




FIG. 15. One loop diagram which contributes to the po- 
larization II = —^GG in the leading order large JV approxi- 
mation. 

The function Z) _1 (s) has a singularity at s — and 
depending on the exact form of N(k) near k — in gen- 
eral, a branch cut starting from the origin in the com- 
plex s plane which may be taken along the negative real 
s axis. Its physical origin is the zero mass Goldstone 
bosons propagating in the internal loop of the polariza- 
tion II in Fig. 15. The discontinuity across this cut 
arises from the non-zero probability for the time depen- 
dent x field to create Goldstone pairs with arbitrarily 
small spatial momentum k above the massless threshold. 
The function D~ 1 (s) also may have one or more zeroes 
in the left-half complex s plane. The zero at 



s = ±4to 2 exp 



32^ 2 



hX B (m 2 



±oo 



(6.17) 



for A# <C 1 is the Landau ghost pole in the far ultraviolet 
which lies outside the range of validity of the large N 
expansion and which in any case does not affect the long- 
time behavior of the inverse Laplace transform. 

In the case of a strictly ther mal d istribution of massless 
bosons ( 3.15 ), the integral in ( |6. 16 ) may be performed in 
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closed form in terms of the digamma function ip and we 
find 



1 



1 

Xr(s 2 /A) ' s 2 



In - 



16tt 2 
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A fl (47r 2 T 2 ) s 



AttT 
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2 



s v \ AttT J 



16tt 2 



2ttT 



4nT 



(6.18) 



A particular case is T = in which case the integral 
vanishes and I? ( ^ 1 (s) possesses a zero at complex s± in 
the left-hand complex plane, 



s± = ±iuj — 7 



(6.19) 



with 



— A p i« 

Co "V 

~ 64^' 



and 

3 

647T 1 



where 



X P i = Xr(s 2 /4 = -ul/4) , 



(6.20) 



(6.21) 



for X p i <C 1. In the frequency ( 6. 20] ) the effective mass of 
the "radial" mode, ignored in the direct quantization of 
the N — 1 massless modes in ( 2.22Q reappears in the large 
TV limit as an oscillation in the real time linear response 
to perturbations about the vacuum (ft = v, x — 0- It may 
be viewed either as this radial degree of freedom of the 
N component $j field, with its effective mass dressed by 
the polarization effects of II in the presence of massless 
Goldstone modes, or alternately and just as correctly, as 
a genuine collective excitation of the composite field \. 
The two descriptions are equivalent since the oscillations 
of y and <f> are constrained by the linearized equations 



(3.7) and (3.£), and there is only one degre e of freedom 
between them. The oscillation frequency (6.20) should 
be compared with the second derivative of the fr ee en - 
ergy potential at T = which vanishes from Eqn. (4.11). 
Thus the "effective" potential is completely ineffective at 
predicting the radial oscillation frequency at zero tem- 
perature. It is clear that a zero temperature the origin of 
the decay rate 7 is the imaginary part of the two-particle 
cut in the polarization diagram illustrated in Fig. 16. 

At finite temperature the analytic structure of -D^ 1 (s) 
in the complex plane with Re s < is as anticipated 
on general physical grounds, i.e. there is a logarithmic 
branch cut beginning at s — along the negative real 
s axis with a discontinuity that monotonically increases 
from zero to its a symptotic value of ±i?i/327r from the 
logarithm in ( 3.15 ) as s runs from to —00. Indeed, the 
digamma function of small argument behaves like 



dlnT(z) 

dz 



1 

IE 

z 



as z -> (6.22) 
is Euler's constant, so that we 



where 7^ = 0.577215 . 
see from the second form of ( 3.18 ) that the logarithms 
cancel and the discontinuity vanishes, leaving a simple 
pole structure, 



% ~ is (6 - 23) 



as s — > 0. Setting this to zero and solving the quadratic 
equat ion informs us that D^, (s) has a zero at s± as in 



( |6.19D with 



J plT — ^r(T' ) 

\r(T 2 ) 



and 



It 



16tt 



T 



(6.24) 



for \s±\ < T and V = 2ire- 2 ^ E T. This condition implies 
that Eqn. (6.24) is only valid in the critical region, 



T c -T 



T 2 



h\ R T? 



(6.25) 



where one would legitimately question the validity of the 
leading order mean field approximation. The approach 
of w 2 ; rp to zero as T — > T c is a manifestation of critical 
slowing down. Notice al so tha t this near critical tempera- 
ture plasma oscillation ( |6.24 ) docs n ot go over smoothly 
to the zero temperature case ( 6.20| ) since opposite lim- 
its of the digamma func tion a re involved, and T — > is 
not permitted in Eqn. (6.24) without violating the as- 
sumption \s±\ <C T. However, it seems likely that the 
simple poles found at zero temperature mi grate contin- 
uously from their locat ion a t s± given by ( 3.20 ) to the 
negative real axis as in (3.24) as the temperature is raised 
from zero to T c . It would be interesting to map out the 
analytic structure of £>y 1 (s) in the left half complex s 
plane as a function of T to study the intermediate tem- 
perature behavior of the linear response function from 
its low temperature form to the critical region, without 
recourse to the asymptotic expansion (6.22). We reit- 



erate that D^ 1 (s) describes the real time oscillations of 
the (f> mean field coupled to the collective plasma oscil- 
lations of the auxiliary field Xj and that the equilibrium 
free energy "effective" potential F has no such real time 
information at either zero or fin ite temperature, as con- 
clusively demonstrated by ( 1.11 ). 

The linear response predictions for the oscillation fre- 
quency and the damping rate in the case of finite tem- 
perature were compared with a numerical evolution with 
Aa = .01 and T = 1. The theoretical predictions of 
uj plT = .09574 and 7t = 1.9894 x 10~ 4 are in very good 
agreement with the numerical results lu p it = .09585 and 
7t = 2x10~ 4 . A similar test was carried out at a larger 
value of the coupling Aa — .1 again with excellent results: 
predictions of u; plT = 0.3028 and 7T = 1.9894 x 10~ 3 
compared with numerical values of u> p it — 0.3031 and 
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7t = 1-9944 x 10~ 3 . The extraction of the behavior of 
the x envelope from the numerical data is described be- 
low. 

Instead of plotting \ directly as a function of time, in- 
formation may be extracted more conveniently by plot- 
ting 



Vt 



\og(y(t + e)) - log(y(t)) 
\og(t + e) - log(t) 



(6.26) 



as a function of time, where each time point t is taken at 
either a peak or a trough of the oscillation and e is the 
time separation between two such neighboring points. At 
late times, if the frequency stabilizes, e goes to a constant. 
Using this fact, alo ng w ith the late time condition t e, 
one can show from ( 3.26 ) that if the % oscillation envelope 
can be fit by an expression of the form 



Ar a cxp -jt 



then at late times, 



Vt 



-7i- 



(6.27) 



(6.28) 



enabling the direct reading off of the power a and the 
exponent 7 as the y-intercept and slope of a straight line 
respectively. It is important to note that this method is 
also a test of whether the frequency is stable as otherwise 
a straight line will not be obtained. The behavior of y t 
against time is shown in Fig. 17 for the two values of the 
coupling constant mentioned above. 
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FIG. 16. The function yt denned in Eqn. (3.26) plotted 
for perturbed thermal initial conditions at two values of Aa- 
The solid lines are the corresponding theoretical predictions. 

When the distribution N(k) is non-thermal the behav- 
ior of Z) _1 (s) is more difficult to analyze, and there is no 
guarantee that it possesses a zero in general. If N(k) is 
peaked at small k then a reas onab le first approximation 
is to neglect the last term in ( 6.16| ). In that case we find 
a zero of D~ 1 (s) at with 



Xpiv 2 
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and 



64tt 



(6.29) 



If the running renormalized Xr is not very small and/or 
N(k) is not sharply peaked at k — then the position 
of this zero of D-\s) will be shifted from jrh29|) , and at 
some values it may even cease to exist. Comparing with 
our numerical data for which kN{k) is shown in Fig. 18, 
we obs erve t hat the approximations leading to the esti- 
mates ( |6.29| ) can only be order of magnitudes at best. 
Indeed the numerical value of the plasma freq uency from 
the data at late times is oj p i = 1.405 whereas (6.2S) gives 
ujpi = v. The damping rate is estimated to be 7 = 0.005u 
(or 0.007?; if the correct uj p i is used), however the enve- 
lope of the oscillations at late times is not even strictly 
exponential as we shall see below. Hence we cannot ex- 
pect to obtain an accurate estimate of the damping rate 
by a simple pole analysis of the Laplace transform D(s). 
Nevertheless, if one fits an exponential to the data at 
late times a value 7 = 0.009u is obtained, agreeing only 
roughly with the estimate. 

Despite the envelope of the oscillations being non- 
exponential, ujpi is very well determined by the regular 
oscillations observed (to a few parts in 10 3 ). To obtain 
a more accurate approximation to the plasma frequency 
consider the following analytic approximation to kN(k), 



kN(k) = 87twq 6 (lu 2 




(6.30) 



where ujq is a free parameter which is close to the numer- 
ical value of the plasma frequency. The accuracy of this 
analytic fit to the data is shown in Fig. 18 for ujq = u) p i. 
Notice that because of dephasing, the N(k) of the initial 
state in the linear response analysis may be identified 
with the N(k) of the particles created earlier in the non- 
linear evolution from the spinodal region (except for very 
small k < 1/t where the dephasing is as yet ineffective) . 
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FIG. 17. Comparison of the the ana lytic approximation 
to the particle number density of Eqn. ( J6.3C| ) (upper curve) 
with the adiabatic number density of the numerical calcula- 
tion (bottom curve). 

Substituting the form ( |S.30 ) into ( |S.16] ) we obtain 
1 v 2 



-D- L (s) = 



Ak(* 2 /4) 



(6.31) 



2uj 2 



dy(l 



,2\ 2 



l^fa 2 
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Because of the small prefactor of the logarithm in 
l/Aij(s 2 /4), the running coupling constant is only a few 
percent different from the bare A near the zero of D^ 1 (s) 
and this running with s can be neglected in lowest order. 
Again we set s = ±iu> p i — 7 and look for a zero of D^ 1 (s) 
near ujq = oA, w ith j/lu <C 1. One could perform the 



integral in ( |6.31 ) exactly but we content ourselves with 
this simple approximation and obtain 



1 v z 
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H — 



dy [(1 - y 2 



(i-v 4 



(6.32) 



Performing the integral we find from the real part of this 
equation, 
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2A 



pi 



A p ; 



(6.33) 



In our simulations where Aa = 1 — \ p i, the output 
value of the calculated plasma oscillation frequency u> p i = 
1.414« is within 0.6% of the measured value lo v \ — 1.405u 
showing the self-consistency of the approximation and 
the fit ( 6.3C ), t o th is accuracy. We also performed the 
integration in (3.16) numerically using the data of Fig. 
18 for kN(k) and reconfirmed w p ; = (1.405±0.001)v from 
the location of the minimum of the real part of D^ 1 (s). 




FIG. 18. The function y t plotted for N(k) = initial 
conditions at two values of Aa- The late time linear fits to 
the data are also shown with the intercepts determining the 
power law prefator. 

The envelope of the x oscillations obtained numerically 
is very well fit by t~ a exp(— jt). This general behavior 
has been checked for different values of the bare cou- 
pling and we have verified tha t it docs not depend on 
which choice of the tOk profiles (C4) and ( |C5| ) one makes 
for the initial conditions. Illustrated in Fig. 19 are the 
Aa = 1 and Aa = .1 cases, where a — 1.5, 7 = .0067 
and a = 1.288, 7 = .00032 respectively. Figs. 19 and 20 



display the complete fit to the damped oscillation along- 
side the actual evolution for the case Aa = 1. Very good 
agreement is evident over long times and large amplitude 
range. The case of Aa = .1 is very similar. 



x 





FIG. 19. The x evoluti on as a function of time (dashed 
line) for Aa = 1. The fit (6.27), shown by the solid curve 
is x(t) = 6.63t -1 ' 5 exp(-.0067t) cos(1.405t + .5). This figure 
shows the fit over an early time range, t = 20, 120. 




FIG. 20. The same \ evolution at late times t = 400, 520 
showing the excellence of the simple fit over long times and 
along with the results of Fig. 19, very good agreement over 
three decades of amplitude. 

As expected from the time dependence of the enve- 
lope of the oscillations observed in Figs. 18, 19, and 20, 
the damping rate cannot be obtained from a simple pole 
analysis. Indeed the imaginary part of D~ 1 (s) evaluated 
from the numerical data does not possess a zero for any 7 
with Dpi fixed at 1.405, although it does decrease mono- 
tonically as 7 is decreased, thus favoring 7 < 0.003 v. 
This value is close to our momentum resolution dk and 
hence we could not go to lower values reliably. In any 
case, the decay of the envelope at late times is due to the 
creation of very low momentum massless bosons which 
have not yet dephased efficiently. Hence it is the behav- 
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ior of the Laplace transform on the branch cut very close 
to s = which determines the late time damping, and 
this behavior depends in turn on the the mode fun ction s 
\fk{t)\ 2 near k = where neither our replacement, ( 5.22 ) 
nor our single pole dominance assumptions are justified. 

In earlier studies of damping a — | power law was en- 
countered and identified as due to the behavior of the 
relevant spectral density function near threshold |Q . In 
the present case there is no simple argument that we are 
aware of which would lead to such a power law at late 
times. 



VII. SUMMARY 

In this paper we have presented a study of nonequilib- 
rium evolution and time dependent behavior of symmetry 
breaking transitions in N component A<3? 4 field theory. 
Starting from an effective action principle for the leading 
order mean field approximation we emphasized the cor- 
respondence of the equations with the Schrodinger evolu- 
tion of a Gaussian density matrix according to a certain 
effective classical Hamiltonian. This is important for em- 
phasizing that no fundamental irreversible behavior has 
been introduced by the mean field approximation. Ex- 
plicit numerical time evolutions from the unstable spin- 
odal region show effective time irreversibility in the form 
of the efficient averaging to zero of the phase informa- 
tion in the density matrix. The effective von Neumann 
entropy of the reduced density matrix measures the infor- 
mation that is lost by discarding this phase information. 
Its general increase with particle creation shows the close 
connection between dephasing and irreversibility. Thus 
our study of effective dissipation and decoherence casts 
some light on fundamental issues in quantum statisti- 
cal mechanics such as the origin of irreversibility, Boltz- 
mann's H-theorem and the quantum to classical transi- 
tion. Here there is much that could be done still in the 
context of the framework presented in this paper, most 
notably to study quantitatively the Poincare recurrence 
cycle(s) expected in a Hamiltonian system and their de- 
pendence on the various parameters of the model. 

In the specific case at hand, <& 4 field theory is a proto- 
type of models of of spontaneous symmetry breaking in 
a wide variety of physical systems. It has been common 
to rely heavily on the thermodynamic free energy and 
equilibrium considerations generally in analyzing these 
systems. One broad and generally applicable conclusion 
of the present study is that this can be very misleading 
where nonequilibrium dynamics is concerned. We have 
shown that the distribution of particles created in the 
mean field time evolution from an initial configuration in 
the spinodal region is generally far from thermal. This 
leads to a final state in which the mean field need not 
be close to minimum of the free energy potential. The 
simple observation that there is a sum rule which the non- 
thcrmally distributed particles can saturate is sufficient 



to resolve this seemingly paradoxical result. 

The linear response analysis of the oscillations about 
the allowed stationary configurations is also a new result 
which is quite different from what consideration of the 
free energy function alone might lead one to expect. We 
have shown that there is a collective plasma mode in the 
radial symmetry breaking direction whose characteristics 
depend on the ambient particle distribution. This also 
should be generally true of the various systems described 
by the same i> 4 field theory. While there is excellent 
agreement in the thermal case for both the plasmon fre- 
quency and damping rate, the non-thermal situation is 
more complicated. In this case, due to the lack of pole 
dominance, the late time behavior cannot be described 
by exponentially damped oscillations. Our numerical re- 
sults for the plasmon oscillation envelope are remarkably 
well-fitted by an expression of the form exp(— jt). 
The plasmon frequency compares well with an analytical 
estimate. However, the precise dependence of the enve- 
lope function on the particle density in the general case 
remains to be more fully investigated. We have outlined 
also how dephasing leads also to effectively classical field 
configurations which one can sample in order to extract 
information on the creation of topological defects during 
the phase transition. This direction is certainly an inter- 
esting one to pursue both in the condensed matter and 
cosmological applications. 

Throughout the paper we have endeavored to bring the 
theoretical framework into close contact with practical 
numerical methods. Indeed one of the main conclusions 
of this work is that the real time dynamics of phase tran- 
sitions can be studied in a concrete way with presently 
available computers. Besides the wide applicability of the 
spontaneously broken $ 4 theory it is interesting for the 
existence of massless Goldstone bosons which are created 
freely during the phase transition. This essential kine- 
matics is shared by systems with an exact gauge symme- 
try such as QCD or general coordinate invariance such 
as gravity. The numerical techniques used in this work 
and the experience gained in treating the massless case 
should prove to be valuable in generalizations to these 
cases. 
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APPENDIX A: THE ENERGY-MOMENTUM 
TENSOR AND ITS RENORMALIZATION 



The effective action (2.7) leads directly to a conserved 



energy-momentum tensor at any order of the large N 
expansion [Q. In this Appendix we give some details 
about the energy-momentum tensor and its renormaliza- 
tion. Because of spatial homogeneity and isotropy the 
only nontrivial components of this tensor are 



(Too) = £ and 



(Al) 



where the energy density e and isotropic pressure p are 
given by 

+ \ J[dk}a k {\f k \ 2 + (k 2 + X )\f k \ 2 } 
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We have used the gap equation ( 2.32 ) in passing to the 
latter expressions in each case. 

The energy-density e contains a quartic but constant 
cut-off dependence {i.e. proportional to A 4 ). This is the 
expected vacuum energy contribution which should be 
subtracted. The conservation of the energy density, 



i = 



(A4) 



is easily checked from the equations of motion, and is not 
affected by the subtraction of the constant quartic diver- 
gence. As it will turn out this single constant subtraction 
is sufficient to yield a cut-off independent conserved en- 
ergy density. 

The cut-off dependence of the isotropic pressure p is 
a bit more complicated and requires a detailed under- 
standing of how Lorentz invariance (or more generally, 
coordinate invariance) is broken by the spatial momen- 
tum cut-off we have introduced in all the mode integra- 
tions. In the first forms of the two expressions (A2) and 
([A3]), i.e. before using the gap equation, there appear 
mode integrals whose cut-off dependences correspond to 
those of a free theory with the mass m 2 replaced by x(t). 
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Divergences in the energy-momentum tensor expectation 
value (Tjlw/) have been studied in the literature by covari- 
ant methods such as dimensional regularization or covari- 
ant point splitting [B5| with the result that these diver- 
gences must be proportional either to the metric tensor 
9nv = diag (— 1, 1, 1, 1) in flat Minkowski spaceti mc o r 
the total derivative "improvement term" of Eqn. (A12) 
below. Because we are using a regulator in the form of a 
cut-off in the spatial momentum integrations which is not 
covariant under Lorentz or general time dependent coor- 
dinate transformations, the quartic and quadratic cut-off 
dependence we shall actually obtain will not have these 
covariant forms. Hence we shall have to perform the 
stress tensor renormalization in a noncovariant manner 
as well to correct for the spurious noncovariant A 4 and A 2 
cut-off dependence, in order to obtain covariant results 
in the end. 

The actual cut-off dependence in both e and p is eas- 
ily analyzed by means of an adiabatic expansion to the 
mode function equation (2.2C ), th e lowest order solution 
of which is given by Eqm (5.12). By substituting this 



lowest order adiabatic approximation to the mode func- 
tions into the mode integrals in the first forms of Eqns. 
(A2) and (A3) we can characterize the most severe de- 



pendence on the ultraviolet cut-off A in the forms, 
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Clearly the A dependence is not such that £o = — Po, as 
required by considerations of general covariance. As has 
been known for some time the reason for this is that the 
spatial momentum cut-off A acts as a noncovariant point 
splitting regulator would, giving terms in the regulated 
(Tfj,u) proportional to S^S^ in which the spatial direc- 
tions i, j = 1, 2, 3 are distinguished. Since such terms do 
not appear in the = v = time component, the en- 
ergy density has the correct A dependence and requires 
no covariantizing correction. However for the pressure 
we should have 



Po 
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on grounds of general covariance. This is easy to enforce 
by hand by adding the difference p' Q — po = —So — po to 
Po above which just corrects for the noncovariant terms 
induced by our momentum cut-off. 



With this prescription to enforce covariance of the 
mode integrals in the pressure, the quartic subtraction re- 
quired on both e and p is the removal of the cosmological 
vacuum energy fi,A 4 /167r 2 , corresponding to a subtrac- 
tive renormalization of the constant cosmological vacuum 
term. We are not allowed to subtract the subleading 
quadra tic and logarithmic divergences of e appearing in 
Eqn. (A5), since, for one thing, they are multiplied by 
the time dependent function x an d hence would spoil en- 
ergy conservation, and for another, such time dependent 
terms do not correspond to a constant vacuum energy. 
Instead we recognize them to be just the correct cut-off 
dependent terms needed to combine with the cut-off de- 
pendent terms in the "classical" energy density, 
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in order to render the total energy density cut-off inde- 
pendent. Because the analogous terms in p are precisely 
the negative of those appearing in e after the correction 
p' a — po has been added to its quantum part, we also 
obtain a (partially) renormalized pressure by the same 
manipulations. 

The reason that we must still perform one additional 
subtraction to obtain a fully renormalized pressure is that 
at the next order i n the WKB adiabatic expansion of 
the mode equation ( [2.20 ), one obtains for the adiabatic 
frequency, 
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If the energy density is calculated to this adiabatic order 
one finds no additional cut-off dependence in e, so no ad- 
ditional subtractions are required for it. Indeed none are 
permitted consistent with the principle of general coordi- 
nate invariance. However, the pressure has an additional 
logarithmic cut-off dependence equal to 
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Since this divergence appears only in the pressure but not 
in the energy it is consistent with the generally covariant 
and conserved form of the "improvement term," pa] 
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in flat space. This term has no effect on the energy den- 
sity for spatially homogeneous mean fields but adds to 
the pressure the total derivative, 
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in which £ is an ar bitra ry parameter. The fact that a 
divergence such as ( |A1C ) appears in the pressure means 
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that we should introduce £ as a free bare parameter of the 
theory from the very beginning, on the same footing as 
the mass p\ and coupling Aa, and allow for the possibility 
that the bare £a 7^ will renormalize in general. Indeed 
this is known to be the case in A$ 4 theory Q and we 
have the renormalization condition, 
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in terms of the j3 function of the coupling constant A. 
We note that the trace of the renormalized energy mo- 
mentum tensor vanishes (for any value of £) in the static 
spontaneously broken vacuum where cf> — x = and 



(A13) 
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where Z x 1 is the same logarithmic renormalization con- This implies that the relativistic equation of state 
s tant as tha t for the coupling constant appearing in Eqns. 
(|2 35|) and fl2.36|) of the text. Thus, 
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and the additional term in the pressure ( A12| ) becomes 
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Comparing the latter expression with (A10) we observe 
that the logarithmic cut-off dependence cancels in the 
sum p + Ap. Hence we must add the improvement term 
and renormalize £ in precisely this way in order to obtain 
a completely cut-off independent pressure. 

To summarize this discussion of the renormalization of 
the energy and pressure we define 
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to be the cut-off independent conserved energy density, 
and 



PR = P - Po - £0 + 
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to be the renormalized pressure of the theory for the 
spatially homogeneous mean fields considered in this pa- 
per. 

Let us also remark that the value £ = 1/6 for the im- 
provement term is also the value chosen for conformal in- 
variance of the scalar theory in the massless limit. Indeed 
the trace of the renormalized energy-momentum tensor 
is 
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h X 2 



v 2 x + 



32tt 2 



1 

Ar 



X 



(A18) 



upon using Eqns. (A5), (A£), and ( A15 ) above. At 
£ = 1/6 the last term vanishes, the first term is the renor- 
malized classical trace of the energy-momentum ten- 
sor and the second term is the one-loop quantum trace 
anomaly, 



holds at late times, independently of the number density 
distribution. The numerically obtained approach to this 
equation of state is shown in Fig. 20. We emphasize 
that this equation of state does not imply relaxation to 
thermal equilibrium since collisional effects are not yet 
taken into account by the leading order large N approxi- 
mation, and the nonthermal nature of the late time state 
is apparent from the created particle distribution of Fig. 
10. 



5=1/6 



FIG. 21. Evolution of the pressure as a function of time 
for £ = 1/6. The approach to the equation of state p — e/3 is 
clearly seen. 



APPENDIX B: 



THE GAUSSIAN DENSITY 
MATRIX p 



In this second Appendix we examine some of the prop- 
erties of the general mixed state Gaussian density matrix 
fl3.16 ) and its time evolution. The discussion will be re- 
stricted mainly to the case of d = spatial dimensions 
to simplify the notation. Generalizations to d > are 
straightforward. First, the time evolution of the density 
matrix is unitary: 



p(t)=U(t)p(0)tf(t), U(t) = 



exp — 1 



(Bl) 
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where H osc is the time-dependent harmonic oscillator 
Hamiltonian. 



H osc (®,V-t) = 



1 



(B2) 



that preserves the Gaussian structure of p under time 
evolution, so that the Liouville equation, 
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is satisfied without taking the trace. 

It is not difficult to find the explicit form of the unitary 
operator U(t) in the coordinate basis, 



(x'\U(t)\x) = (2irihv(t))~i x 
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in terms of the two linearly independent solutions to the 
classical evolution equation, 



dt 2 



x 2 (t) 
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u(0) = v(0) = 1 



(B5) 



u(0) = v(0) = . 

This same U(t) also evolves the quantum operators, 

*(t) = U\t) $(Q) U(t) = 4>(t) + af(t) + t/»(t) 

P(t) = U\t) P(0) U(t) = p(t) + af(t) + t/*(f) . (B6) 

Mathematically, the three Fock space bilinear opera- 
tors aa, and a^a + aa^ generate the Lie algebra of 
the symplectic group Sp(2) = SU(1, 1) ^ SL(2, R) which 
is the group of homogeneous linear transformations of 
phase space (<&, P) which preserves the anti-symmetric 
classical Poisson bracket {#,P} = 1, i.e. 



P' 



a b 
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Since this is one condition on four real parameters the 
group Sp(2) is a three parameter group. If the a and 
Heisenberg operators are appended to these three, the 
algebra again closes upon itself, forming a five parame- 
ter group, the inhomogencous mc ta plcct ic group, IMp{2) 
PH . The unita ry ev olution (Bl), (B4) of the Gaussian 
density matrix (3.16) is an explicit representation of this 
group's action. 

The form of the density matrix in the time- independent 
Heisenberg basis is quite easy to obtain from the form of 
the the transition matrix clement (x\n) given by Eqn. 



(5.7) of the text. By substituting the integral represen- 
tation of the Hermite polynomials, 
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where C is a closed contour around the origin of the com- 
plex z plane into the expression 
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and interchanging the orders of integration we can per- 
form the double Gaussian integral over the shifted vector 
(x — 0, x' — <t>) first. Using the standard formula 
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the exponent of the resulting expression simplifies con- 
siderably. Letting z = e 10 and z'* = e~* e we are left 
with 
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Expanding the last exponent in a Taylor series we find 
that only the terms with n' — n survive with Eqn. (5.S) 
the final result. As discussed earlier in Section VI, in this 
Heisenberg basis the density matrix is time-i ndepe ndent 
and diagonal. This fact allows the writing of ( B14 ) in an 
operator form, 
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as far as computing matrix elements in this Heisenberg 
basis is concerned. 

By making a different IMp{2) group transformation it 
is also possible to diagonalize ( |B2[ ) at any given time, 
bringing the quadratic Hamiltonian into the standard 
harmonic oscillator form, H osc — 4^ (aa^ + a^a) with 
a time dependent. This adiabatic particle basis is re- 
lated to the time-independent Heisenberg basis by the 
relations, 



f = af + Pf* , 
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with / the adiabatic mode function defined by Eqn. 
( 5.12 ) of the text, and 



a=y*{j-iuf) , 
[3=-~f{f + iujf) , 

-4>f) ■ 

The Bogoliubov coefficients a and (3 obey 



1 



(B17) 



(B18) 



and therefore may be expressed in terms of a real param- 
eter 7 and two phases, 



a = cosh 7 e 1 ^ , 
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The density matrix in the adiabatic particle number ba- 
sis can be expressed in terms of these time dependent 
parameters at any given time t. The direct evaluation of 
the expression analogous to ( B10[ ) in the h number basis 
is quite tedious, and is accomplished most rapidly by use 
of a coherent state basis as discussed in Ref. [^9| . Making 
the identifications, 
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and x ~~ * 7 m the notation of Ref. p9|| , Eqn. (5.32) of 
that work yields the desired matrix element in the case 
of zero mean field, namely, 
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where P™ is an associated Legendre polynomial. 

The first important feature of this expression for our 
purposes is the phase factor when n ^ n' . If the exact 
mode function / is rewritten in the form 
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this phase 
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Thus, even in the adiabatic limit, where ^- < 1, the 
phase angle 9 depends linearly on ti me a nd the off- 
diagonal elements of the density matrix ( B22| ) are rapidly 
varying i n tim e. On the other hand, the diagonal ele- 
ments of (B22) are independent of this phase angle and 
consequently much more slowly varying functions of time. 
Indeed, in the case of zero mean field d> = 0, the adiabatic 



invariant W of Eqn. ( 5.20 ) is 

W = N - N = a sinh 2 7 (B25) 
and (n|p|n) depends only upon a (a constant) and 7: 
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Now if cr > 1 the Legendre polynomial is not necessarily 
positive which means that we cannot interpret the diago- 
nal matrix elements of the density matrix in the adiabatic 
particle number basis as a classical probability distribu- 
tion in the general mixed state case. However, if cr = 1 
then the argument of the Legendre polynomial vanishes 
for any 7^0. Since 



P2t{0) = (-1) 



■ (21- 1)!! 



(B27) 



for n = 21 even but P 2 e+i (0) = for h = 2t + 1 odd, the 
diagonal matrix element (B26) simplifies considerably in 
the pure state case, 



(n = 2£\p\n = 21)1 



(2^-1)!! , l2£ 
sedry tanh 7 



2 e e\ 



with the mean number of created particles, 
N = Smh 7 = 2hu ' 



(B28) 



(B29) 



which are Eqns. ( 5.25 ) and ( 5.26 ) of the text. These 
results were reported in ||. 

In the pure state case the even n diagonal matrix ele- 
ments of the density matrix are positive definite and may 
be interpreted as the probabilities for observing £ uncor- 
related particle pairs in the adiabatic particle basis. The 
odd n diagonal matrix elements vanish since particles can 
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only be created in pairs from the vacuum. Otherwise (i.e. 
for a > f), the muc h mo re complicated and non-positive 
definite expression (B26) shows that there is no simple 
classical probability interpretation for the density ma- 
trix. The restriction to zero mean fields </> = </> = is 
not serious for spa tially homogeneous backgrounds be- 
cause of Eqn. (3.31) which shows that all the mean fields 
vanish except for the k = mode. Hence, at late times 
all the k > modes may be treated classically if dephas- 
ing is effective and typical classical field configurations in 
t he en semble can be constructed as in Eqns. (5.30) and 
(5.31 ) of the text, provided the k = mode is excluded. 



APPENDIX C: NUMERICAL METHODS 



In order to solve the system of equations ( 2.15 ), ( 2.20 ), 
and ( p4] ) as an initial value problem, we have to spec- 
ify the initial conditions of the mean field <$>, its time 
derivative and the mode function and its time derivative. 



Then we have to solve the gap equation ( 2.34 ) at the 
initial time tg ■ In order to have a finite set of renormal- 
ized equations we have to choose the mode functions so 
that the high momentum modes coincide with the zeroth 
order adiabatic vacuum described by 



A (to) 



A (to) 



1 



y/2u)k(to) ' 



2u> k (t ) 



/fe(to) 



(CI) 
(C2) 



with u>l(to) = k 2 + x(to)- It is easy to verify that for 
initial conditions with (f>(to) — 0, it is also true that 
x(to) = by inspecting the time derivative of the gap 
equation. This simplifies the form of /fc(to) to 



A (to) = -iu k {t )fk(t ) 



(C3) 



For an initial state with positive square effective mass 
X we could use the same form of the mode functions 
of (CI) also for the low momentum modes. However, 
if we wish to investigate the case of a "quench" within 
the unstable spinodal region by initial conditions with a 
negative x(0), we have to modify the initial values of the 
low momentum modes in order to avoid the singularity 
of /fe(to) = l/V^fc at k 2 = — x(to). We have used two 
different profiles of the frequency Uk for the initial mode 
functions: 



^(to) = fc 2 +x(to)tanhf^±fM 

V lx(*o)| 

u;Ut ) = k 2 + X (to)exp(-v 4 /k 4 ) . 



(C4) 
(C5) 



At large momentum these profiles coincide with the adi- 
abatic vacuum frequency as required. In all cases, except 
for thermal initial conditions, the initial number of quasi- 
particles, N(k) = 0. 



Numerical simulations were performed on a massively 
parallel computer using a momentum grid with 32000 
modes. The upper cut-off was set at A = 5v, implying a 
grid resolution: 



dk 



A 



32000 



The field </> was scaled in units of v so that k is also 
given in units of v , x m units of v 2 , and the time t in 
units of v^ 1 . Most of the results discussed here are for a 
bare coupling constant Aa = 1 but other values were also 
investigated. The corresponding renormalized coupling 
con stant s at the renormalization point |x(to)| a- re given 



by (2.35) with m replaced by |x(to)|- With Aa = 1, and 



for <j>(to)/v = 0, \ R = 0.990036, while for </)(t )/v = 0.5, 
A7? = 0.990248. The energy densities in these cases are 
e = 0.1237 and e = 0.06978 respectively. 

The mode equations were stepped forward in time us- 
ing a sixth-order adaptive time-step Runge-Kutta inte- 
grator. The time-steps were controlled by tracking the 
evolution of x- Energy conservation to parts per million 
was achieved over the temporal range of typical evolu- 
tions. 
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